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ABSTRACT 


The  object  of  this  thesis  is  to  develop  a  procedure 
for  calculating  the  dynamic  displacement  response  of 
the  platform  of  a  fixed  offshore  structure  acted  upon  by 
a  regular  wave  train.  The  structure  considered  has  four 
legs  in  a  square  configuration,  with  the  waves  impinging 
normal  to  one  side  of  the  square „  The  procedure  may  be 
manipulated  for  use  with  other  leg  configurations  and  wave 
directions . 

The  thesis  is  written  in  a  manner  useful  to  a  designer 
of  fixed  offshore  structures .  An  iterative  procedure  may 
be  used  to  arrive  at  the  critical  wave -displacement 
combination. 

The  types  of  waves  considered  are  those  for  which  Stokes • 
3rd  approximation  applies.  A  modified  version  of  Morison's 
theory  is  used  for  drag  forces.  Families  of  curves  are  pre- 
sented for  use  in  predicting  these  forces  and  their  centroids 
The  curves,  plotted  in  dimensionless  form,  show  the  limits 
of  applicability  of  Stokes8  3rd  approximation.  MacCamy- 
Fuchs  diffraction  theory  is  used  for  inertia  forces. 

To  calculate  the  dynamic  response  of  the  platform  the 
authors  fit  the  vibration  problem  of  the  structure  to  the 
classical  theory  for  a  linear,  single -degree -of -freedom 
system.  The  wave  forces  are  expressed  in  terms  of  a  Fourier 
series.  Displacements  are  calculated  for  twelve  positions 
in  one  wave  cycle,  and  a  displacement  curve  is  drawn.  A 
simple  tabular  form  is  presented  for  calculating  these  dis- 
placements. 

Computer  programs  (in  FORTRAN)  are  given  for  performing 
the  bulk  of  the  force  and  displacement  calculations  men- 
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tioned  above.  The  theoretical  procedure  for  unbraced  and 
braced  leg  configurations  is  tested  by  experiment  for  32 
cases.  For  these  tests  the  ratio  of  maximum  theoretical 
to  experimental  displacements  varies  from  about  0.5  to  1.7o 

It  is  concluded  that  the  designer  should  incorporate 
a  factor  of  ignorance  of  3  in  his  design. 

The  authors  show  that  the  largest  waves  are  not 
necessarily  the  most  critical  for  the  structure.  Rather, 
the  critical  design  wave  is  a  function  of  the  frequency 
ratio  (wave  frequency/structure  natural  frequency).  The 
designer  must  consider  a  range  of  significant  waves  which 
give  frequency  ratios  of  about  1.1  or  less. 

It  is  recommended  that  calculations  be  performed  for 
some  real  structures  in  ocean  waves.  From  a  comparison 
of  the  results  with  observed  data  a  revised  factor  of  ig- 
norance might  be  obtained. 

In  Appendix  I  an  approach  to  the  solution  for  vortex 
shedding  forces  and  frequencies  in  a  train  of  surface 
waves  is  presented.  Existing  vortex  theories  are  modified 
to  apply  to  a  vertical  cylinder  in  finite  amplitude  waves, 
and  boundaries  are  set  on  the  predicted  maximum  forces . 
Experiments  are  conducted  from  which  vortex  forces  are  seen 
to  fall  within  the  boundary  limits  set  by  theory.  The  ex- 
periments also  validate  the  proposed  theory  that  the  shedding 
frequency  is  a  predictable  integral  multiple  of  wave  fre- 
quency. It  is  recommended  that  a  procedure  be  developed  for 
determining  platform  displacements  due  to  these  forces. 


Thesis  Supervisors   Donald  R.F.  Harleman 

Title:  Associate  Professor  of  Hydraulic! 
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I.   INTRODUCTION 

The  problem  of  wave-induced  vibrations  in  fixed 
offshore  structures  suddenly  was  cast  into  the  public 
news  focus  when  on  January  15,  196l,  the  United  States 
Air  Force  Texas  Tower  Number  Four  collapsed  80  miles 
off  the  New  Jersey  coast  in  a  heavy  storm  with  the  at- 

r  -\* 

tendant  loss  of  28  lives  1.1] . 

Prior  to  that  incident,  and  since,  no  one  has 
computed  the  dynamic  displacements  of  such  a  structure 
acted  upon  by  wave  forces,  to  the  best  of  the  authors' 
knowledge.  The  design  approach  has  been  to  make  the 
natural  frequency  of  the  structure  relatively  high 
(about  50  cycles  per  minute)  and  then  simply  treat  the 
wave  force  as  if  it  were  statically  applied.  This 
method  is  crude  at  best,  and  can  be  very  misleading. 

During  the  winter  of  1958-59  a  platform  motion 
study  [2]  was  made  of  Texas  Tower  No.  4.  It  was  ob~ 
served  that  the  tower  platform  displacements  were 
greater  for  10  and  11  ft.  waves  (about  ±J>   in.)  than 
for  ^O  ft.  waves.   Obviously  the  highest  wave  was  not 
the  most  critical  for  this  structure. 


Numbers  in  brackets  refer  to  literature  citations, 
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The  Air  Force  was  interested  in  platform  motions 
of  Texas  Tower  No.  4  because  they  wanted  to  know  what 
errors  might  be  introduced  into  the  installed  aircraft - 
tracking  radar  system,   (For  this  particular  situation 
it  was  felt  that  rotational  motion  would  probably  be 
more  of  a  problem  than  translational  motion.) 

This  thesis  develops  a  method  of  calculating  the 
dynamic  translation  displacements  of  a  fixed  offshore 
structure  in  a  regular  wave  train.  The  theory  developed 
has  been  experimentally  checked  by  the  authors.  The 
experimental  model  was  patterned  somewhat  similar  to  the 
United  States  Coast  Guard's  Buzzards  Bay  Light  Station, 
located  off  the  coast  of  Massachusetts  f^] . 

The  Coast  Guard  has  initiated  a  long  range  pro- 
gram of  replacing  lightships  by  fixed  offshore  structures. 
The  two  completed  to  date  (Buzzards  Bay  and  Brenton  Reef 
Light  Stations)  are  braced  structures  with  four  legs. 
The  U.S.  Navy  built  a  similar  but  considerably  larger 
structure  at  Argus  Island,  near  Bermuda.  The  U.S.  Air 
Force  built  three  Texas  Towers,  which  were  three-legged 
structures.  Many  other  fixed  offshore  structures  have 
been  built  in  the  past  also,  particularly  in  connection 
with  drilling  and  mining  operations  in  the  Gulf  of  Mexico. 
The  potential  exists  [4  J  for  many  more  offshore  oil- 
drilling  rigs  to  be  built  around  the  world  in  the  years 
to  come. 


There  has  been  at  least  one  case  reported  J_5  J 
in  which  vortex  shedding  was  important  in  waves.  In 
that  case  after  a  one  week  storm  a  two  foot  diameter 
pile  suffered  a  fatigue  failure  caused  by  transverse 
vibrations  of  2.5  seconds  period  in  12  feet  high 
waves  of  13  seconds  period.  The  authors  of  this  thesis 
suspected  that  vortex  shedding  forces  might  have  a 
noticeable  effect  on  offshore  structures.  Some  investi- 
gation of  this  subject  has  been  included  in  this  thesis. 
The  authors  concluded  from  their  literature  survey  that 
no  previous  work  has  been  done  on  the  question  of  vortex 
shedding  forces  on  a  stationary  vertical  cylinder  in  a 
wave  train. 


II.   PROCEDURE 

A,  Theoretical  Procedure 
1.  Hydrodynamics 
General 

This  section  deals  only  with  those  forces  on  a 
cylindrical  pile  acting  in  a  direction  parallel  to  that 
of  wave  propagation.  These  forces  are  termed  longitudinal 
forces  to  distinguish  them  from  the  transverse  or  "lift" 
forces  caused  by  vortex  shedding.  The  transverse  forces, 
although  requiring  consideration  in  a  practical  structure 
design,  caused  such  small  displacements  of  the  model  under 
study  that  it  was  not  feasible  to  obtain  experimental 
verification  of  an  analytic  approach.  Because  of  this, 
the  transverse  force  investigation,  both  theoretical  and 
experimental,  is  included  as  Appendix  I  and  is  not  dis- 
cussed further  in  this  section. 

The  ultimate  design  procedure  would  be  to  make  use 
of  a  statistical  analysis  of  ocean  wave  and  energy 
spectra  in  the  development  of  the  input  forcing  functions 
for  the  structure  displacement  problem. 


Steps  have  been  taken  in  this  direction,  for  example 
[6],   although  the  problems  encountered,  not  only  in  ob- 
taining a  sufficient  scope  of  reliable  data  [ 7  J,  but  also 
in  applying  this  data  to  the  development  of  a  useful  de- 
sign criteria,  have  yet  to  be  solved.  For  this  reason 
the  authors  are  restricting  the  investigation  to  a 
regular  train  of  surface  waves  whose  forces  are  capable 
of  relatively  simple  mathematical  superposition. 

In  the  past  a  typical  design  procedure  has  been  to 
decide  upon  a  "design  wave,"  which  is  the  largest  wave 
to  be  expected  (within  specified  confidence  limits )  over 
the  planned  lifetime  of  the  structure.  Using  applicable 
wave  theory,  one  may  then  determine  the  wave  forces  and 
moments  acting  over  an  incremental  length.  These  may  be 
summed  up  in  a  tabular  form  to  give  the  total  force  and 
moment.  [ 8,  3  J.  Some  designs,  also,  have  considered 
a  breaking  wave  of  size  somewhat  smaller  than  the  design 
wave  [9,  10 J. 

If  the  dynamics  of  the  problem  are  considered  it 
is  obvious  that  this  design  wave  may  not  be  the  wave 
which  causes  the  greatest  excursions  (or  stresses)  in 
the  structure.  Rather  it  is  to  be  expected  that  a  wave 
(of  significent  height)  whose  fundamental  or  some  low 
order  harmonic  corresponds  closely  to  the  natural  fre- 
quency of  the  structure  may  be  the  critical  wave  from 
a  design  standpoint.  This  possibility  has  been  Indicated 


where,  over  an  extended  period  of  observation,  a  7.5  foot 
wave  caused  greater  platform  motion  of  Texas  Tower  No.  2 
than  much  larger  waves  did./~9l    A  similar  observation 
is  cited  in  the  Introduction  for  Texas  Tower  No.  4. 

A  serious  situation  also  may  exist  when  the  wave 
length  is  approximately  equal  to  the  leg  spacing,  for 
in  this  case  the  fundamental  wave  forces  (and  possibly 
some  harmonics)  will  be  in  phase. 

The  proposed  procedure,  therefore,  presents  a 
method  whereby  several  waves  of  various  lengths  and 
steepness  may  be  considered  over  the  range  of  possible 
water  depths.  The  critical  displacement -wave  combination 
is  then  arrived  at  by  an  iterative  process.  To  allow  for 
a  solution  by  this  method  within  a  reasonable  length  of 
time  families  of  curves  are  developed  such  that  total 
force  and  moment  may  be  determined  by  a  simple  multi- 
plication. This  is  intended  to  supplant  the  tedious  in- 
cremental tabulation  and  summation  mentioned  previously. 
The  curves  (Figures  III,  IV,  V)  are  presented  in  terms 
of  dimensionless  parameters  such  that  they  may  be  applied 
to  any  situation  for  which  the  proposed  theory  is  appli- 
cable. 

A  pile  spacing  greater  than  10  pile  diameters  has 
been  found  experimentally  to  be  large  enough  such  that 
proximity  effects  are  negligible.  /ll7    This  spacing 
is  exceeded,  not  only  in  the  vast  majority  of  offshore 
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structures  of  the  type  considered,  but  also  for  the 
model  which  was  the  object  of  the  authors'  experi- 
mental study. 

Wave  Theory 

H  »  wave  height,  measured  from 
crest  to  trough. 

L  *  wave  length. 

d  -  water  depth. 

m     ■  height  of  wave  crest  above 
still  water  level. 

Several  theories  on  wave  motions  have  been  developed 
and  tested  experimentally  with  varying  degrees  of 
success.  These  theories  are  found  to  be  applicable 
over  finite  ranges  of  wave  parameters  and  water  depths, 
with  no  one  theory  satisfying  all  conditions.  Since 
the  waves  of  interest  in  this  study  are  of  the  re- 
latively shallow  water  finite  amplitude  type,  two 
theories  seem  to  be  of  particular  interest  -  namely, 
Stokes'  third  approximation  and  the  Solitary  wave  theory 
These  theories  are  discussed  in  detail  elsewhere  (for 
example  [4,  12,  15,  14"J)  and  will  not  be  repeated  here. 
However,  for  purposes  of  evaluating  this  study,  the 
following  limits  of  applicability  may  be  noted  from 
these  references. 


Stokes1  3rd  approximation: 

a.)   C^q/H)  ^  0.625 
b.)   (H/L)   ^  0.142 

Solitary  wave: 

a. )  waves  at  or  near  breaking  point 
b.)   (d/L)  — >•  0 

Theoretical  limits  on  ocean  wave  parameters: 

a.)  (%/H)  ^  0.75 

b.)  (H/L)    ^  0.142 

c.)  Hmax  -  O.78   (d-(H-^Q))  for  shallow  water 

d.)  from  the  above,  (H/d)  ^  O.653 

Using  the  maximum  value  of  (/7?  /H)  from  Stokes1 
theory  and  the  value  of  H    above,  one  sees  that 
(H/d)    is  O.603  by  Stokes'  theory.  There  is  a 
gap  between  this  value  and  O.653  for  which  the  soli- 
tary wave  theory  would  appear  to  be  more  applicable. 
However,  waves  have  been  observed  breaking  on  Martha's 
Vineyard  with  values  of  H/d  *  O.588.  [15] 

Because  this  present  investigation  must  be  limited, 
it  was  decided  to  use  only  Stokes'  3rd  approximation.  For 
the  solitary  wave  and  other  theories  the  method  of  approach 
in  determining  forces  and  force  centroids  basically  would 
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be  the  same,  with  the  curve  families  of  this  thesis  en- 
larged to  include  their  particular  ranges.  Bretschneider 
£l4j  has  developed  graphs  which  appear  to  be  a  good  start- 
ing point  for  determining  the  range  of  applicability  of 
the  several  theories,  as  well  as  values  of  H,  L,  Oh   and 
maximum  crest  velocity  (Um^)  for  design  waves.  These 
values  could  serve  as  first  trial  inputs  for  the  pro- 
cedure to  be  presented  here.  For  a  conservative  estimate 
of  the  force  exerted  by  a  breaking  wave  on  a  piling  a 
value  of  4  times  the  force  as  determined  from  Stokes' 
theory  for  a  wave  of  the  same  height  may  be  used.  /"l6l 

For  reference  purposes  the  applicable  formulas 
from  Stokes"  3rd  approximation  are  included  here  in 
dimensionless  form.  f_ljj  The  reader  is  referred  to 

the  list  of  symbols  at  the  beginning  of  this  thesis  and 
Figure  I  for  use  with  these  equations. 

Let  C  •  wave  celerity,  ft. /sec. 

J  -  2ird/L  (1) 

The  vertical  displacement  of  a  particle  from 
its  mean  position  is  given  by  (2) 


d 


*u  fraij  yd)J sine  -  (i  u  'gn^J cos  2j 


(2) 


direction  of 
wave  travel 


Still  water  level 


i    * 

<: P 1  X  Y 


}'yf 


if      Bottorn_ 


Figure  I.  Definition  Sketch  of  a  Wave 


At  the  wave  crest  (2)  becomes 


max  -   1   3  slnh(J  S/d)   3  ir   sinh(2J  S/d ) 


(3) 


and  at  the  trough * 


y 


trough   _ 
3~^  "  2 


1   f  -slnh(J  S/d)   g  tt  slnh(2J  S/d)  c 
d7H  (    sinh  J      H  ITH   (si^  J}4 


At  the  still  water  level  (S/d  -  l),  (2)  gives  the  wave 
profile 
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From  (5),  with  0  -  90°  the  height  of  the  crest  above  the 
still  water  level  is 

and  for  9  ■  270°  the  distance  to  the  trough  is 

%'E         1    f   ,  .  3  ir   sinh(2J)   5  ,7x 

The  wave  celerity  is 

tE^-Vt***  jf1  +   (T7g)2I2(00Bh  2J)2+  2CSSh  2J^7?    (Q) 
y-^-       |/J  (  *Z7H'  8(sinh  J)4  JJ 

The  particle  velocity  at  any  depth  S/d  is 
At  the  wave  crest    (9)  becomes 
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u 
max 


I*    ir  _C_  C  cosh  (J  S/d)    t     1     /  U  cosh(2J  S/d)  N 
■=  ^,:rair?   slnhJ     ^'(sinh  J)2("2^  <.lnh  J)2  V 


k/gd  'JgcT 


♦wissw?  ™ 


(slnh  J) 
At  the  wave  trough  (9)  becomes 


^trough    t  _C__  C  cosh  (J  S/d)    t     1     /  13  cosh  (2  J  S/d)x 

+  *  1  j-cosh(2JS^d)j>  (n) 

L/H  2  L  (sinh  J)2  -'^ 


Force  Theory 

Using  Stokes1  3rd  approximation  for  wave  motions, 
D.R.F.  Harleman  and  W.C.  Shapiro  of  the  Massachusetts  In- 
stitute of  Technology  Hydrodynamics  Laboratory  have  de- 
veloped a  force  theory  which  shows  good  agreement  with  ex- 
periment . £  17 1  In  developing  this  theory  they  applied 
certain  modifications  to  Morison's  Force  Theory  for  drag 
and  the  MacCamy-Fuchs  Diffraction  Theory  for  inertia  forces. 
In  this  thesis  the  authors  will  utilize  the  approach  as  pre- 
sented by  Harleman  and  Shapiro.  Families  of  curves  will  be 
developed  from  their  equations,  these  curves  to  be  used  for 
the   iterative   design  procedure  mentioned  previously. 
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As  will  be  shown  the  only  forces  of  interest  are 
the  maximum  values  of  the  drag  and  inertia  components  for 
both  crest  and  trough  regions. 

FD    «  maximum  drag  force  in  crest  region 
(FD),  «  maximum  drag  force  in  trough  region 

FT  ■  maximum  inertia  force  in  crest  region 
(FT ),  -  maximum  inertia  force  in  trough  region 
D  »   cylinder  diameter 

(PD)t  ,-rDd2CD  trough  B  (13) 


l+7</d 


where 


A~  J      ~w^-  d(s'/d)  <lZl> 

o 
B   i=        f  U   tgU6ho   d(S'/d)  (15) 


and  S1  is  the  vertical  distance  to  the  actual  particle 
position  and  is  therefore  equal  to  the  mean  particle 
position  S  plus  or  minus  the  vertical  displacement  y. 
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For  the  crest  S'  -  S  +  y  (16) 

max 

and  for  the  trough      S»  -  S  -  ytrough  (17) 

It  can  be  seen  that  (14)  and  (15)  do  not  lend  themselves 
to  regular  integration,  but  must  be  solved  graphically  or 
by  some  other  means.  The  steady  state  drag  coefficients 
CD  and  CD  4-r0Uo.h  to  be  used  in  (12)  and  (1J>) ,   respectively, 
are  determined  as  functions  of  Reynolds  numbers,  defined 
as  follows: 


(18) 


„,   y  u  max  p.   D  \l    g  d  A~ 

-       -  t trough   .  D  |   S  5 ■  B  (19) 

trough      ^        ^  W     %-H 


where  \/  u     and  v/  u  .    .  are  the  root  mean  squares 
Y        max     y    trough  H 

of  the  velocity  distributions  for  the  crest  and  trough, 
respectively. 


j 


u 


2    _    /  A  g  d2 
max     V  d  + 


% 


(20) 


V 


u2       -    /  B  g  d2  ,-x 

trough    v/  d  +  (%j-H)  *   ' 

In  (18)  and  (19),  A  and  B  are  determined  from  (14) 
and  (15),  respectively,  and  %/d,  -*§ —  from  (6)  and  (7), 
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Figure  11 ,  Steady  state  drag  coefficient  versus 

Reynolds  number  for  circular  cylinders, 
fl7,  18J. 

respectively.  Figure  II  is  a  plot  of  steady  state  CD  vs 
for  circular  cylinders. 


FT  -  it  ^  WJ)  t  H  sinli  I  cos  9 
I   4  cosh  J  , 


(22) 


2 

(Ft)*  *  )/   Tv.Pt  "  Slnh  Z    C0S  9 

1  I't   4  cosh  J 


where  I  -  J  +  j~75   sin  e 


(23) 
(24) 
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and  9  is  that  angle  of  wave  position  (0°  <  9  <-  180° 
for  (22),  and  180  5  9  ^  360  for  (25) )  such  that  P-.  and 
(Fj)t  are  maximized,  respectively.  Because  9  will  be 
very  nearly  equal  to  0°  and  180°  for  these  situations, 
it  can  be  seen  from  (24)  that  little  error  (and  great 
simplification)  is  introduced  by  letting  9  «  0°  and 
l30°  in  (22)  and  (23 ),  respectively. 

This  gives    Fj  -  -(Fx)t  -  — ^  %D — -  tanh  J      (25) 


Let  S1  -  center  of  action  of  the  force  Fp  and 
(S> ).  -  center  of  action  of  the  force  (F-)t 


1 


+  7o/d 


Then 


|1.      1      s,/d  ^dd(s./d)  (26) 


1+i£ 


■H 


-^i-  r  S'/d^gg**2  d(S'/d)  (27) 


ST  -  center*  of  action  of  force  F-j- 

5I  -  /gfcffh  i  [  1+J81I*J-  cosh  JJ        (28) 
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Harleman  and  Shapiro  have  shown  that  S'  can  be 
assumed  constant  for  0  ^~  0  -  180°,  and  (S')t  can  be 
assumed  constant  for  180°  -  9  *  j560°  with  little  loss  in 
accuracy.  Also,  with  the  assumption  that  S  ,  is  equal 
for  crest  and  trough  by  (28),  it  is  reasonable  to  assume 
that  it  is  constant  for  any  value  of  0. 

3y  use  of  (5,7,10,11,14,15,26,27)  values  of  A  and 
B  for  (12,15,18,19)  and  |^-  and  ^^  for  (26)  and  (27) 
were  computed  for  families  of  the  parameters  H/L  and  d/L 
on  an  IBM  7090  digital  computer. 

The  results  are  plotted  in  Figures  III,  IV,  and  V 
in  dimensionless  form.  Values  of  H/d  equal  to  0%6Q3, 
which  was  mentioned  as  the  approximate  limit  of  applica- 
bility of  Stokes1  theory,  and  0.652,  which  is  the  maximum 
value  this  parameter  may  obtain  (page  8  ),  are  also 
shown  in  these  figures.  These  values  indicate  the  boundary 
for  use  of  the  figures,  as  well  as  the  boundary  up  to  which 
the  solitary  wave  theory  could  be  used  to  extend  the  range 
of  applicability ,  The  computer  program,  representative 

output,  and  comments  are  included  as  Appendix  D.  In  ad- 
s'1      s"f 
dition  to  A,  B,  g-  ,  and  (g-)t>  the  following  are  given 


as  outputs: 


% 
d 

d 
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H/d  =-  0.653 

1/d  «  0..6QI 


0.1 


603  (limit 
Stokes ■  theory ) 

U/d  <*   0.653  (theoretical 
limit) 


Drag  Force  Centroids  for  Wave  Height/Wave  Length 
(H/L)  and  Water  Depth/Wave  Length  (d/L)  Ratios. 

(S  /d)  -  centroid/depth  ratio  for  creat 

(S  trougft  /d)  "  centroid/depth  ratio  for  trough 
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Design  Force  Procedure 

It  has  been  shown/"  17]  that  drag  force  is  a 

2 
function  of  sin  9,  while,  for  the  assumption  that  F- 

occurs  at  9  »  0°,  180°,  inertia  force  is  a  function 
of  cos  9.  The  forcing  functions  for  use  in  the  dif- 
ferential equations  of  platform  displacements  will  there- 
fore be  of  the  following  form: 

0°  £  9  -   180°  2 

F  -  FD  sin  9  +  F,-  cos  9        (29) 

180°  ^  9  *  560°  B  _  /T?  x  flin2 


F  =  (FD)t  sin*9  +  F  cos  9       (50) 


For  convenience  the  formulas  and  procedures  used 
in  the  determination  of  these  forces  and  their  points  of 
application  are  summarized  below. 
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Pd  ,      t    Y2  °D  A  (12) 

(PD)t  -  ^Dd2C2D  trouSh  B  (15) 

Fl  «    **  f  H     tanhJ  (25) 


Cn  from  Figure  II  for    /R  »  ?  V      s  d  A  (18) 

v    '  1  + 


/r-    ix'um  rigure    ix    lor     n\  ■  — 


Vd 


S  trough  from  F1*ure  II  **    *  -f  y/     g%°-H         <ld' 

A  from  Figure  IV 
B  from  Figure  V 


a~    2 


w*  ♦**«$>? 


(sinh  J) 
or  from  computer  output  (not  plotted) 


»H 


V-V-i  <»> 


si    and    (s8)t    from  Flgure  in 


5i  '  rra  j  / *  +  JsinhJ- cosh  JJ      <28> 
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2.  Vibrations 


General 


Once  the  wave  forces  on  the  legs  have  been  deter- 
mined the  next  problem  is  to  calculate  the  resulting 
displacement  of  the  platform  as  a  function  of  time,  or 
wave  position  9  with  respect  to  the  front  (upstream) 
legs  of  the  structure. 

This  thesis  considers  only  the  longitudinal  dis- 
placement x  of  the  platform,  taken  as  positive  in  the 
direction  of  wave  propagation.  It  has  been  observed 
that  the  transverse  vibrations  in  the  experimental  runs 
were  relatively  negligible,  with  an  amplitude  of  perhaps 
three  percent  of  that  of  the  longitudinal  vibrations : 
For  the  dimensions  of  real  structures  the  transverse 
vibrations  may  be  as  important  as  the  longitudinal 
vibrations. 

This  thesis  does  not  consider  wind  force  and  its 
effect  on  platform  displacement.  Certainly  with  large 
waves  there  will  also  be  high  wind,  and  the  wind-caused 
displacement  must  be  added  to  that  induced  by  the  waves. 
In  a  steady  wind  the  problem  is  one  of  statics ,  not 
vibrations,  unless  vortex  shedding  is  important.  A  useful 
discussion  of  wind  loads  is  given  in  [4J  .  Water  currents 
can  be  treated  in  the  same  way. 
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The  structure  considered  has  a  square  platform  of 
weight  W,  supported  on  four  (N-4)  legs  of  length  X  and 
diameter  D,  arranged  in  a  square  configuration  of  leg 
spacing  h  on  a  side.  The  case  of  no  leg  bracing  is  con- 
sidered first,  after  which  the  case  of  added  bracing  is 
considered,, 

The  investigation  is  limited  to  the  case  of  the 
waves  impinging  normal  to  one  face  of  a  square  con- 
figuration of  legs.  Once  the  reader  understands  the 
concept  of  this  procedure  he  should  be  able  to  manipulate 
the  procedure  to  fit  any  practical  leg  configuration  and 
a  limited  number  of  angles  of  wave  impingement .  When  the 
direction  of  wave  propagation  is  not  parallel  to  an  axis 
of  symmetry  of  the  braced  leg  configuration,  the  calcu- 
lations become  complicated  by  the  facts  that  generally 
the  platform  translation  is  not  in  the  direction  of  the 
applied  force  and  that  rotation  of  the  platform  may  be- 
come significant . 

The  legs  of  the  thesis  model  are  fixed  at  the  base. 
This  condition  was  assumed  for  the  Argus  Island  structure 
ML 9l  .  In  the  case  of  Texas  Tower  No,  4  the  footings 
beneath  each  leg  were  25  feet  in  diameter  filled  with 
cement,  and  sunk  or  embedded  into  the  ocean  floor  to  a 
depth  of  18  feet .  After  the  collapse  divers  reported 
that  the  footings  were  apparently  in  good  condition  with- 
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out  any  evidence  of  fracture,  movement,  or  scour  [20J  . 
The  opinion  was  that  the  footings  were  not  cocked  out  of 
a  perpendicular  plane  with  the  ocean  floor  [21J  .    In 
other  words,  the  assumption  of  fixed  footings  was  valid 
for  Texas  Tower  No.  4. 

The  legs  of  the  thesis  model  are  fixed  at  the  plat- 
form,, This  assumption  was  made  in  the  design  of  Texas 
Tower  No.  4c  Later  investigation  [2J  indicated  that  the 
actual  restraint  condition  was  closer  to  being  pin-ended 
than  fixed. 

As  a  practical  matter  in  the  consideration  of 
vibrations,  it  should  be  pointed  out  that  large  quantities 
of  surging  oil  or  water  in  tanks  may  amplify  tower  motion. 
This  problem  generally  can  be  reduced  by  installing 
several  small  tanks  in  lieu  of  one  large  tank  and  by  in- 
stalling Internal  baffling  in  the  tanks . 

Reference  ([20J  on  the  collapse  of  Texas  Tower  No.  4 
provides  much  useful  and  interesting  food -for -thought 
for  the  designer  of  fixed  offshore  structures. 

The  reader  interested  in  experiments  with  models 
should  realize  that  it  is  not  possible  to  scale  model 
results  to  fit  the  prototype.  The  problem  is  that  drag 
and  inertia  components  of  the  wave  force  do  not  scale 
in  the  same  manner.  Consequently  the  reader  should 
center  his  interest  in  the  concept  of  this  thesis,  and 
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not  on  the  numerical  results  for  any  particular 
run. 

Classical  Vibration  Theory 

The  approach  considered  will  be  to  fit  the 
vibration  problem  of  the  structure  to  the  classical 
theory  for  a  linear,  single -degree -of -freedom  system. 
There  are  many  sources  in  the  literature  where  one 
may  find  this  theory.  The  authors  have  chosen  to 
follow  the  notation  of  Den  Hartog  [221. 

From  Newton's  Law  one  can  arrive  at  the  follow- 
ing differential  equation  of  forces: 

mx  +  Gx  +  kx  «  P  sin  u)t  (52) 

where   x  »  displacement,  a  function  of  time 
x  *  dx/dt  =  velocity 
x*  =  crx/dtr  *  acceleration 
m  -  mass  with  single  degree  of  freedom 
G  ■  damping  coefficient 
k  ■  spring  constant 
p  *  magnitude  of  applied  force 
co  *=  circular  frequency  of  applied  force 
t  »  time 

The  steady  state  solution  to  (32)  can  be  written  as 

x  .  xQ  sin(60t  -  0)  (33) 

where  x  and  0   are  as  yet  unknown,  and : 
o  ,      * 
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x  =  displacement  amplitude 
0  *  phase  angle  by  which  x  lags  P  sin  cot 


Insert  (33)  into  (32)  to  obtain: 


x~ ?0  (34) 


/(CO))2  +  (k  -  mco2) 


2  ' 


0  m   arc  tan  p  (35) 

From  consideration  of  free  vibration  without  damp 
ing  one  can  find  that  the  undamped  natural  frequency 
Wn  is: 


^n"/*  (56) 


From  consideration  of  free  vibration  with  damping 

one  can  define  critical  damping  C  as: 

c 


Cc  =  2  Jink     «  2mCOn  (37) 


For  values  of  damping  C  less  than  C  the  free 
vibration  is  truly  oscillatory.  For  values  of  C 
greater  than  C  the  free  "vibration"  is  not  oscillatory, 
but  rather  the  mass  creeps  back  exponentially  to  its 
equilibrium  position  from  an  initial  displacement. 
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Insert  (36)  and  (37)  into  (3*0  and  (35)  to  obtain: 


'STATIC 


n  j 


(38) 


where  x 


STATIC 


(39) 


JZf  -  arc  tan 


V°n 

l-(f)2 
n 


(40) 


It  will  prove  useful  to  have  numerical  tables  for 

expressions  (38)  and  (40)  for  entering  arguments  of  ~ 

n 
and  small  values  of  C/c  .  Table  I  is  presented  for 


values  of 


Static 


Table  II  is  presented  for  values  of  0, 


Unbraced  Structure 

Mass  m  must  represent  the  mass  of  the  platform 
and  a  certain  portion  of  the  mass  of  the  legs.  By 
applying  Rayleigh's  Energy  Method  to  this  case  (Appendix 
A)  it  is  possible  to  compute: 


m 


l(w  +M«) 

gv    35   ' 


(41) 
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TABLE  I  -  MAGNIFICATION  FACTORS  (58) 


TABLE  II  -  PHASE  ANGLES  (degrees)  (40) 


">n  X 

0 

0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.25 

1.07 

1.07 

1.07 

1.07 

1.07 

1.07 

1.07 

1.07 

1.07 

0.50 

1.55 

1.53 

1.33 

1.33 

1.33 

1.33 

1.33 

1.33 

1.33 

0.70 

1.96 

1.96 

1.96 

1.96 

1.95 

1.94 

1.94 

1.93 

1.92 

0.80 

2.78 

2.78 

2.77 

2.76 

2.74 

2.72 

2.68 

2.65 

2.62 

0.85 

5.60 

5.59 

3.58 

3.55 

3.50 

3.45 

3.39 

3.31 

3.24 

0.90 

5.26 

5.24 

5.19 

5.06 

4.92 

4.76 

4.58 

4.39 

4.20 

0.95 

10.21 

10.05 

9.54 

8.85 

8.09 

7.35 

6.69 

6.07 

5.54 

1.00 

00 

50 

25 

16.7 

12.5 

10.0 

8.33 

7.14 

6.25 

1.05 

9.81 

9.61 

9.09 

8.35 

7.59 

6.84 

6.18 

5.60 

5.09 

1.10 

4.76 

4.74 

4.66 

4.55 

4.40 

4.22 

4.04 

3.84 

3.66 

1.15 

5.10 

3.09 

3.07 

3.04 

2.98 

2.92 

2.84 

2.77 

2.69 

1.20 

2.27 

2.27 

2.26 

2.24 

2.22 

2.19 

2.16 

2.12 

2.08 

1.50 

1.45 

1.45 

1.45 

1.44 

1.43 

1.43 

1.42 

1.40 

1.39 

1.50 

0.80 

0.80 

0.80 

0.80 

0.80 

0.80 

0.79 

0.79 

0.79 

1.75 

0.48 

0.48 

0.48 

0.48 

0.48 

0.48 

0.48 

0.48 

0.48 

2.00 

0.55 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

5.00 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

5.00 

0.04 

0.04 

o.o4 

0.04 

0.04 

0.04 

0.04 

0.04 

0.04 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.25 

0 

0 

1 

1 

1 

2 

2 

2 

2 

0.50 

0 

1 

2 

2 

3 

4 

5 

5 

6 

0.70 

0 

2 

3 

5 

6 

8 

9 

11 

12 

0.80 

0 

3 

5 

8 

10 

12 

15 

17 

20 

O.85 

0 

4 

7 

10 

14 

17 

20 

23 

26 

0.90 

0 

5 

11 

16 

21 

25 

30 

34 

37 

0.95 

0 

11 

21 

30 

38 

44 

49 

54 

57 

1.00 

0/0 

90 

90 

90 

90 

90 

90 

90 

90 

1.05 

180 

168 

158 

148 

140 

134 

129 

125 

121 

1.10 

180 

174 

168 

162 

157 

152 

148 

144 

140 

1.15 

180 

176 

172 

168 

164 

160 

157 

153 

150 

1.20 

180 

177 

174 

171 

168 

165 

162 

159 

156 

1.30 

180 

178 

176 

174 

171 

169 

167 

165 

163 

1.50 

180 

179 

177 

176 

174 

173 

172 

170 

I69 

1.75 

180 

179 

178 

177 

176 

175 

174 

173 

172 

2.00 

180 

179 

178 

178 

177 

176 

175 

175 

174 

3.00 

180 

180 

179 

179 

178 

178 

177 

177 

177 

5.00 

180 

180 

180 

179 

179 

179 

179 

178 

178 
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where  g  =  gravitational  constant,  normally 
52.2  ft/sec. 2 
w  ■  weight/unit  length  of  leg 
N  =  number  of  legs 


If  the  configuration  had  been  such  that  the  legs 
were  considered  to  be  pinned  to  the  platform  the 
fraction  to  be  used  in  the  second  term  of  (4l)  would 
be  y$/\h§    (Appendix  A). 

It  is  shown  in  Appendix  A  that  the  added  mass 
effect  (of  accelerated  water  mass)  is  negligible  in 
still  water .  For  forced  vibrations  there  is  the  further 
consideration  that  the  water  particles  are  moving  with 
wave  orbital  velocities „  For  this  latter  situation  it 
is  simply  assumed  that  added  mass  is  negligible .  Later 
in  the  experimental  results  (Appendix  G)  it  will  be 
seen  that  this  assumption  is  valid. 

For  damping  C  it  is  necessary  to  resort  to  experi- 

Q 

ment  to  determine  values  of  •=-   for  the  plastic  model  of 

Cc 
this  thesis.  The  damping  consists  of  structural  damping 

in  the  plastic  and  viscous  damping  by  the  water.   Again 

from  Den  Hartog  [22  If or  the  case  of  free  vibrations  with 

viscous  damping  one  finds  that  the  logarithmic  decrement 

d         fits  the  following  relationships : 

of  a  loge  ^  (42) 
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where  x.,  and  x  are  two  successive  positive  (or 
negative)  decaying  displacement  amplitudes. 

/  -      g    ,.f-  (^3) 

i  -  (§  ) 

cc 

For  values  of  C/C  up  to  8#,  and  somewhat  larger, 
the  following  approximation  is  "exact"  within  sliderule 
accuracy: 

§-  *  SF^  «*> 

With  this  expression  one  can  calculate  values  of 
structural  damping  by  observing  free  vibrations  in  air, 
assuming  that  the  air  causes  negligible  damping. 

For  completeness  at  this  time  it  is  noted  that: 


co  -      q  (45) 

' i1 

1  -  (§-) 
°c 

where   q  •  damped  natural  frequency 

With  this  expression  one  can  calculate  experimental 

values  of  00  . 
n 

It  is  shown  in  Appendix  A  that  in  still  water  for 
these  particular  experiments  the  theoretical  viscous 
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damping  is  11$  maximum  of  the  structural  damping.   In 
the  steady  state  vibrations  there  is  the  further  con- 
sideration that  the  velocity  associated  with  viscous 
damping  is  not  x  but  is  the  velocity  relative  to  the 
water  particles,  which  are  moving  with  wave  orbital 
velocities.  For  these  reasons  viscous  damping  by  the 
water  is  neglected.  The  experimental  results  (Appendix 
F)  support  this  action. 

In  Appendix  A  the  following  expression  is  derived 
for  spring  constant  k: 


(46) 


12NEI 

where 

E  ■  Young's  modulus 

and 

I  *  moment  of  inert 

with  bending  of  one  leg 

Actually  k  must  be  determined  experimentally  for 
the  plastic  model  of  this  thesis  because  a  unique 
value  of  E  cannot  be  obtained  from  a  handbook  for  the 
plastic. 

From  (36,41,46); 


CO  .      12gN?Lw0  (47) 
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The  forcing  function  in  (J2)  is  P  sinOJt. 
The  expression  for  the  wave  force  is  not  this  simple, 
but  if  it  can  be  represented  by  a  series  of  terms 
of  the  form  P  sin  a)  t  one  can  find  an  equal  number 
of  terms  for  x,  each  with  its  own  phase  angle.  By 
superposition  of  these  terras  one  can  calculate  the 
platform  displacement  x. 

The  forces  acting  on  each  leg  are  equal  in  magni- 
tude. The  forces  on  the  front  legs  are  in  phase 
(because  of  the  direction  of  wave  impingement  on  the 
structure).  The  forces  on  the  rear  legs  are  likewise 
in  phade  with  each  other,  but  they  lag  the  forces  on 
the  front  legs  by  a  lag  angle  ^RL,  which  one  can  deduce 
to  be: 

_/     36Q  h  degrees   2rrh  radians       /kox 
^RL     L  "  L  l   ' 

For  the  drag  force  acting  on  one  leg  it  has  been 

stated  (29,30)  that  in  the  range  0°  £=.   0  ^  180° 

p 
the  force  can  be  expressed  as  FD  sin  co  t.  At  9  «  90° 

the  line  of  action  of  this  force  is  at  a  height  s1  above 
the  bottom.  For  180°  4=.    9  ^  3600  the  drag  force  is 
(FD^trouKh  sin  ^  *'  which  acts  at  a  height  (s'L.  The 
inertia  force  acting  on  one  leg  can  be  expressed  as  FT 
cos  co  t,  acting  at  a  height  s,  when  9  =  0°,  180°. 
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Only  a  fraction  of  these  forces  can  be  considered 
to  act  directly  on  the  platform.  Let  P  represent 
the  force  required  at  the  level  of  the  platform  to 
produce  the  same  static  deflection  as  force  P  acting 

at  height  a  on  a  leg  of  length  X   •  An  influence 

p  _ 

fraction  o/p  has  been  computed  in  Appendix  A,  with 

the  numerical  results  given  in  Figure  VI. 

It  has  been  stated  earlier  that  FD  can  be  con- 
sidered to  act  at  I"'  throughout  the  range  0°  ^=  9  ^  180° 
This  means  that  in  this  range; 


a   s1 

— -  a  — 

Use  Figure  VI,  and  write 

P 


t 


(F  is  a  positive  number.) 
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(*9) 


F-l  -  FD.^  (50) 

With  similar  reasoning; 


(51) 


p2  ■  -(fdVp?  (52) 


FIGURE  VI  -  INFLUENCE  FRACTION 

~p — r — 

Fraction  ~  of  Force  P  to  be  Considered  Acting  On 
Platform  When  P  Acts  at  Height  a  on  an  Unbraced 

Leg  of  Length  jt. 
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7 "  f  (55) 

F^  «  Pj.pS  (54) 


Now  one  can  write  expressions  for  the  force  acting 
on  the  platform  due  to  the  force  acting  on  one  front 
leg. 

force  -  F.,  sin2cot  +  F,  cos  00  t  (0%u>t  4=.   180°)   (55) 
force  «  -F2sin260t  +  F,  cos  lo  t   (l80°4  uJt  4z  360°)   (56) 

.  A  Fourier  series  is  calculated  in  Appendix  A  to  re- 
present this  force.  For  the  cases  considered  in  this 
thesis  only  five  terras  are  significant.  The  following 
expression  results; 


F(t)  -  1  g  2  +  —(F-L  +  F2)  sin  cot 

-  -jr- (F.,  +  F2 )  sin^  a>  t  +  F,  cos  oo  t 

"  *  X  4  2^  cos2  60t  (57) 


The  angle  ^L  has  been  defined  (48)  for  a  harmonic 
force  varying  as  cot,  One  can  deduce  that  if  the  force 
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varies  as  n  60  t  the  appropriate  lag  angle  is 

lag  angle  =  n0RL  (58) 

At  this  point  all  of  the  necessary  theory  is  avail 
able  for  computing  displacement  x  as  a  function  of  9 
at  the  front  legs .   One  has  simply  to  set  up  a  con- 
venient method  for  a  hand  numerical  calculation. 
Fortunately  a  simple  tabular  method  can  be  found.  For 
purposes  of  illustration  it  seems  best  to  show  this 
method  in  conjunction  with  a  sample  calculation.  The 
reader  is  therefore  referred  to  Appendix  B. 

The  numerical  calculation  can  also  be  performed 
on  a  digital  computer.  Such  a  program  was  used  in  the 
preparation  of  the  Results,  Section  III.  The  program 
is  given  in  Appendix  D. 

Braced  Structure 

The  solution  to  the  problem  for  a  braced 
structure  may  be  approached  by  applying  rough  modifi- 
cations to  the  solution  for  the  unbraced  structure. 

For  mass  m  it  was  found  previously  (41 )  by 
Rayleigh's  Energy  Method  that  13/35  of  the  mass  of  the 
legs  should  be  added  to  the  mass  of  the  platform. 
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Assume  that  the  same  ratio  holds  for  the  bracing  and 
legs  of  the  braced  structure.  Define  the  ratio  r^  as 


m   mass  of  legs  +  mass  of  bracing         /,-q} 
m  *     mass  of  legs  *   ' 


rn 


The  assumption  is : 


m  -  i(W   +  i|  rmN  wi>  )  (60) 


From   (36,46,60) ; 


a,     -    /   -     ^i (61) 

In  the  design  of  the  Argus  Island  structure  the 
fraction  ~  was  used  in  lieu  of  ^.  j_  23J    It  appears 
that  this  fraction  was  assumed,  not  derived. 

Damping  for  the  plastic  model  must  be  determined 
experimentally,  as  before  (44) . 

r,In  most  metal  structures,  nearly 
all  the  damping  observed  can  be 

•  attributed  to  the  technique  of 
fabrication.  For  example,  a  good 
weld  can  be  distinguished  from  a 
poor  one  by  its  effect  on  the  part's 
ringing  after  impact:  the  poorer  the 
joint,  the  more  damping.  Riveted 
joints  generally  contribute  more 
damping  than  welded  ones,  and  bolted 
joints  more  than  riveted  ones . H  * 

Hamme,  R.N.,  [24,  p.  37-2J 
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"Bare  steel  structures,  particularly 
welded  rather  than  riveted,  have  very 
low  damping  ranging  from  perhaps  1.0 
to  3.0  per  cent  of  critical  damping."  * 


If  a  designer  is  in  doubt  as  to  what  damping  exists 
in  his  structure  he  will  generally  be  well  advised  to 
assume  the  smallest  reasonable  value  in  order  to  be  on 
the  "safe  side." 

Stiffness  k  remains  as  expressed  in  (46)  for  the  thesis 
model  because  no  bracing  was  added  in  planes  where  it  would 
have  stiffened  the  structure  for  resistance  to  longitudinal 
forces.  If  such  bracing  had  been  added,  experimental  dis- 
placements would  not  have  been  large  enough  to  be  ob- 
servable o 

For  a  normal  braced  structure  one  must  compute  ap- 
propriate values  for  stiffness  k  and  influence  fraction 
p 
p^.   A  good  textbook  on  this  matter  is  [_25J  .  Probably 

P 
two  values  of  _o  will  be  sufficient  to  fair  a  curve,  be- 

P 

cause  the  curve  will  be  very  similar  in  shape  to  that  of 

Figure  VI. 

The  case  of  Texas  Tower  No.  4  provides  an  interesting 

lesson  in  engineering  with  regard  to  stiffness  furnished 

by  the  bracing.  On  this  tower  the  bracing  consisted  of 

pin  connected  tension  members.  It  was  concluded  [_2~J 

that  the  underwater  bracing  was  essentially  ineffective 


Housner,  G.W.,f  24,  p.  50-2Q~? 
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for  platform  excursions  up  to  +3  inches  because  of  the 
large  clearances  in  the  pin  connections.  The  lesson  is 
that  joints  should  be  welded,  which  is  the  method  usually 
used. 

It  is  seen  in  (12,13)  that  drag  force  varies  directly 

with  leg  diameter  D.  It  is  seen  in  (25)  that  inertia 

p 
force  varies  directly  as  D  .  For  each  of  these  forces 

the  line  of  action  is  independent  of  D.  Assume  that  if 

bracing  existed  whose  long  dimension  was  parallel  to 

the  direction  of  wave  propagation  the  wave  forces  on 

this  bracing  would  be  negligible  compared  to  the  forces 

on  the  bracing  whose  long  dimension  was  parallel  to  the 
wave  crests .  Let  J^  be  the  cumulative  length  of  all  of 

a 

the  hydrodynamically  significant  underwater  bracing 
(below  the  still  water  level),  of  diameter  DB«  For 
simplicity  assume  that  the  wave  forces  on  a  length 
d(d  ■  water  depth)  of  bracing  are  the  same  as  if  the 
bracing  were  uniformly  distributed  in  one  additional  leg 
of  diameter  D^.  These  assumptions  lead  to  the  following 
relationships  for  forces  acting  on  the  braced  structure: 

Vb  "     W     <F1>B  "     VF1>    <Vb  "     VF2  <62) 

Where  r»-i+^^)  .(63) 
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(FZ)B-  rrFr  (P5)B-  r^  (64) 

where     r  ~  1  +  ^T^)  (65) 


With  these  modified  values  for  the  forces  to  be 
used  in  (57)  one  can  make  a  hand  calculation  for  dis- 
placement x  as  in  Appendix  B}  or  a  computer  calculation 
as  in  Appendix  D. 

Vibration  Stresses 

Up  to  this  point  no  mention  has  been  made  of 
stresses .  It  should  be  obvious  that  if  platform  dis- 
placement is  greater  for  the  dynamic  situation  than  for 
the  static  situation  the  attendant  vibration  stresses 
throughout  the  structure  will  also  be  greater  by  the 
same  ratio .   (Deflections  are  "small/1  and  the  stress- 
strain  relationship  is  linear  elastic.) 

The  maximum  dynamic  displacement  has  been  ob- 
tained from  the  plots  of  x  vs  9  that  resulted  from 
carrying  out  the  foregoing  procedure.  As  a  result  of 
the  calculation  for  stiffness  k  there  is  available  a 
moment  (and  force)  distribution  (which  is  simply  re- 
lated to  a  stress  distribution  by  simple  beam  theory). 
The  stress  distribution  is  for  a  selected  force,  which 
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results  in  a  calculated  platform  displacement.  This 
static  displacement  is  the  denominator  of  the  sought - 
for  ratio.  All  values  in  the  static  stress  distri- 
bution are  multiplied  by  this  ratio  to  arrive  at  the 
vibration  stress  distribution. 

In  addition  to  the  vibration  stresses  one  must 
also  take  into  account  static  stresses  due  to  such 
causes  as  wind  forces,  current  forces,  column  loads, 
thermal  effects,  fabrication,  etc.  These  stresses 
can  be  superposed  to  arrive  at  the  total  stress  at  any 
point . 

Bo  Experimental  Procedure 

The  experimental  work  of  this  thesis  was  carried 
out  in  the  Wave  Tank  of  the  Hydrodynamics  Laboratory 
of  Massachusetts  Institute  of  Technology  in  March  1962. 

The  primary  objective  of  the  experiments  was  to 
measure  platform  displacement  of  a  model  of  a  fixed  off- 
shore structure  under  the  influence  of  a  regular  train 
of  waves.  Thirty -two  runs  were  made.  In  addition,  runs 
were  made  to  measure  the  longitudinal  and  transverse 
wave  forces  on  a  single  vertical  cylinder.  More  will  be 
said  later  about  the  various  runs  and  the  data  that  was 
collected. 
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Most  of  the  equipment  used,  with  the  exception  of 
the  models,  was  in  existence  prior  to  these  particular 
experiments  and  has  been  described  in  detail  in 
several  references  J  17,26,27,28  J  . 

Wave  Tank  and  Wave  Generator 

Figure  VII  is  a  photograph  of  the  wave  tank, 
which  is  90  feet  long,  2j-  feet  wide,  and  3  feet  deep. 
Basically  it  is  a  steel  structure,  but  the  walls  and 
40  feet  of  the  bottom  are  plate  glass. 

At  the  near  end  of  the  tank  there  is  a  hinged 
plate  wave  generator,  hinged  to  the  tank  bottom.  A 
rod  connects  the  upper  end  of  the  plate  to  an  eceentrlc 
point  on  a  rotating  wheel.  Wave  height  is  governed  by 
adjusting  the  eccentricity.  Wave  period  is  governed 
by  adjusting  the  speed  of  rotation  of  the  wheel.  The 
speed  can  be  determined  by  a  photocell  arrangement 
wired  to  an  electronic  counter. 

At  the  far  end  of  the  tank  there  is  a  25  foot 
long  beach  arrangement,  whose  purpose  is  to  absorb  energy 
and  prevent  (or  reduce)  undesirable  wave  reflections. 
For  these  experiments,  however,  it  was  decided  to  avoid 
wave  reflection  altogether  by  starting  the  wave  generator 
for  each  run  and  recording  data  before  the  reflected  wave 
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Figure  VII  -  Photograph  of  Wave  Tank 


Figure  VIII  -  Photograph  of  Braced  Structure  Model 
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had  time  to  return  to  the  model.  After  some  trial  and 
error  with  the  waves  to  be  used  it  was  decided  that  data 
would  be  recorded  between  the  passing  of  the  eighth 
and  ninth  wave  crests  at  the  test  stand  for  all  runs. 

Test  Stand 

The  test  stand  is  seen  in  Figure  VII.  It  is 
located  approximately  40  feet  downstream  from  the  wave 
generator.  It  consists  of  a  triangular  frame  of  steel 
angles  supported  by  three  legs  of  5-inch  pipe  straddling 
the  tank.  The  stand  is  guyed  to  the  floor  by  steel 
cables  to  increase  its  rigidity.  A  steel  angle  was 
added  to  connect  the  stand  to  the  laboratory  wall  to 
improve  transverse  rigidity,  with  a  consequent  improve- 
ment in  the  rigidity  previously  supplied  by  the  cables. 

Models 

Figure  VIII  is  a  photograph  of  the  model  of  a 
braced  structure.  It  was  fabricated  of  plexiglas  by 
Forest  Products,  Inc.,  of  Cambridge,  Mass. 

The  legs  are  J  inch  (actually  0.505  inch  average) 
rod  42  inches  long,  with  16  inches  on  a  side  between 
leg  centers.  The  leg  diameter  was  decided  upon  as  a 
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compromise  between  the  requirments  for  a  large  diameter 
for  adequate  experimental  forces  and  for  a  small  diameter 
for  adequate  flexibility  in  order  to  observe  deflection. 
The  bracing  is  £  inch  rod,  with  horizontal  members 
located  14  and  28  inches  above  the  base.  The  platform 
is  5/8  inch  thick  (actually  0.610  inch  average)  by 
l8-l/l6  inches  square.  A  5/8  inch  bolt  in  the  center 
of  the  platform  serves  as  a  pointer  for  displacement 
measurements  and  holds  the  weights  (to  be  discussed  later) 
The  base  of  the  model  is  J  inch  thick  by  20-1/8  inches 
by  22  inches,  with  beveled  leading  and  trailing  edges. 
The  model  also  was  used  in  the  unbraced  condition. 

Single  Cylinders 

In  order  to  measure  longitudinal  and  transverse 
wave  forces  two  different  aluminum  cylinders  were  used. 
They  were  attached  to  the  portal  gage  by  a  flange  ar- 
rangement, as  seen  in  the  photograph  of  Figure  IX.  The 
i  inch  cylinder  cleared  the  tank  bottom  by  1/8  inch. 
The  1  inch  cylinder  cleared  the  tank  bottom  by  12-5/8 
inches. 
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Figure  IX  -  Photograph  of  Single  Cylinder  in  Position 
for  Measuring  Transverse  Forces 
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Movie  Equipment 

A  16  mm  Bell  and  Howell  movie  camera  was  used 
to  record  platform  motion  and  time,  as  seen  in  the 
photographs  of  Figures  X  and  XI.  Camera  speeds  of 
24  and  48  frames/second  were  used. 

Recorder 

A  four-channel  direct  writing  model  150  San- 
born recorder  was  used  to  record  wave  profiles  and 
force  traces.  A  paper  speed  of  25  mm/sec  was  found 
to  be  most  satisfactory.  Time  marks  could  be  made 
in  the  paper  margin  upon  signal  from  a  manual  push- 
button. 

Wave  Profile  Measurement 

Wave  profiles  were  measured  by  a  resistance 
wave  gage,  which  can  be  seen  in  Figures  IX  and  X. 
The  gage  consists  of  two  platinum  wires  (0.008  inch 
diameter)  insulated  from  each  other  and  oriented 
vertically  with  a  spacing  between  them  of  J  inch. 
When  the  wires  are  partially  immersed  in  water  an  al- 
ternating current  flows  which  is  essentially  directly 
proportional  to  the  depth  of  immersion.  The  current 
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Figure  X  -  Photograph  of  Braced  Structure  Model 
in  Position  for  Measuring  Displacement 


Figure  XI  -  Photograph  Showing  Scheme  for  Recording 
Platform  Displacement 
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produces  a  signal  which  is  amplified  and  recorded  by 
the  Sanborn  recorder. 

The  gage  can  be  raised  and  lowered  by  a  graduated 
rack  and  pinion  arrangement.  It  is  calibrated  statically 
by  raising  and  lowering  it  through  known  distances  in 
still  water. 

Two  gages  are  used  to  determine  wave  length.  The 
gages  are  raised  until  their  probes  are  immersed  about 
J  inch  in  the  wave  crest .  The  distance  between  the 
gages  is  adjusted  by  trial  and  error  to  one  wave  length, 
at  which  time  two  crest  marks  will  be  made  simultaneously 
on  the  Sanborn  recorder. 

Force  Measurement 

Horizontal  forces  were  measured  by  a  portal 
gage,  seen  in  Figure  IX  directly  above  the  cylinder 
flange  and  oriented  in  this  case  to  measure  transverse 
forces .  The  gage  measures  horizontal  force  in  terms  of 
shear  deflection.  Its  design  is  such  as  to  make  the  out- 
put essentially  independent  of  the  distance  to  the  point 
of  application  of  the  force,  or  bending  moment.   A 
Schaevitz  linear  variable  differential  transformer  is 
used  to  convert  the  shear  deflection  of  the  bottom  plate 
of  the  gage  into  an  electrical  signal,  which  is  amplified 
and  recorded  by  the  Sanborn  recorder. 
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The  force  gage  can  be  calibrated  in  the  positive 
and  negative  directions  by  the  use  of  weights  suspended 
on  a  string  which  is  run  over  a  pulley.  Two  pulleys 
are  seen  in  position  on  the  sides  of  the  tank  in 
Figure  IX. 

Model  Runs 

Water  depth  was  maintained  at  27  inches  through- 
out the  experiments. 

Four  different  waves  were  used,  with  the  charac- 
teristics as  listed  in  Table  III. 


TABLE  III 

Experimental  Wave  Characteristics 

Wave 

H(feet) 

L(feet) 

d/L 

H/L 

A 

0.240 

15.23 

0.148 

0.0158 

B 

0.313 

10.86 

0.207 

0.0289 

C 

O.361 

8.82 

0.255 

0.0409 

D 

0.404 

7.33 

0.307 

0.0552 

The  model  was  used  in  two  configurations,  un- 
braced and  braced.  For  each  configuration  and  wave  com- 
bination four  structural  natural  frequencies  were  in- 
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vestigated  by  adding  weights  to  the  platform.   The 
weights  used  were  0,3,6,  and  10  pounds. 

A  total  of  32  model  runs  was  made.   For  each  run 
platform  displacement  and  time  were  recorded  by  the 
movie  camera,  and  wave  profile  was  recorded  by  the 
Sanborn  recorder.   A  time  mark  was  made  oh  the  re- 
corder paper  at  a  known  time  as  read  on  the  clock  that 
appeared  in  the  movie. 

Single  Cylinder  Runs 

Single  cylinder  runs  were  made  for  each  wave 
on  a  J  inch  cylinder  to  measure  longitudinal  and  trans- 
verse forces.  Runs  were  made  for  each  wave  on  a  1  inch 
cylinder  to  obtain  additional  data  on  transverse  forces, 

For  each  run  wave  profile  and  force  trace  were 
recorded. 

Miscellaneous  Measurements 

Movies  were  taken  of  the  free  oscillations  of 
the  unbraced  and  braced  structures  in  air  and  in  water 
for  all  combinations  of  the  four  added  weights.   Dis- 
placement and  time  were  recorded  for  three  cycles. 
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An  experimental  value  of  spring  constant  k  was 
determined. 

Experimental  values  of  the  influence  fraction 

p 

p-°-  (Figure  VI)  were  determined  as  a  function  of  j. 

The  weight  of  the  plastic  in  the  model  was  re- 
corded. 

Water  temperature  was  recorded. 
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III.  RESULTS 

Figure  XII  presents  the  theoretical  and  experi- 
mental results  for  platform  longitudinal  displacement  x 

The  terminology  "Run  AO"  indicates  a  run  with 
wave  A  and  0  lb.  of  added  weight  on  the  platform  of  the 
unbraced  model.  The  letter  "Z"  specifies  the  braced 
model. 

Experimental  wave  characteristics  are  given  in 
Table  III. 

c 

Values  of  structural  damping  jr-  used  in  the  cal- 

c 
cul at ions  are  given  in  Table  X. 
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Figure  Xlla  -  Platform  Displacements  for  Wave  A 
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Figure  Xllb  -  Platform  Displacements  for  Wave  A 


0.20 
Platform 
Displ . 
in 
Inches 


-0.20 


0.20 


-0.20 


Theory  

Experiment  000 


Run  A3  Unbraced 


90  180  270 

Wave  Phase  Angle  at  Front  Legs  in  Degrees 


Run  A3Z  Braced 


0.626 


560 


VCH  4/19/62 


56 


;:""■*.>  .".:*. "rr"".'r::T*  ".r"'"r '.".'";'"  —-;-••— 


Figure  XIIc  -  Platform  Displacements  for  Wave  A 
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Figure  Xlld  -  Platform  Displacements  for  Wave  A 
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Figure  Xlle  -  Platform  Displacements  for  Wave  B 
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Figure  Xllf  -  Platform  Displacements  for  Wave  B 
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Figure  Xllg  -  Platform  Displacements  for  Wave  B 
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Figure  Xllh  -  Platform  Displacements  for  Wave  B 
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Figure  Xlli  -  Platform  Displacements  for  Wave  C 
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Figure  XII j  -  Platform  Displacements  for  Wave  C 
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Figure  Xllk  -  Platform  Displacements  for  Wave  C 
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Figure  XIII  -  Platform  Displacements  for  Wave  C 


0,20 

Platform 
Displ. 
in 
Inches 


-0,20 


Theory  

Experiment  ooo 


Run  CIO  Unbraced 


1.452 


0,40 


90  180  270 

Wave  Phase  Angle  at  Front  Legs  in  Degrees 


0,20 


0 


■0.20 


0 

0 

0 

Run  C10Z  Braced 
-£--  1.416 


360 


VCH  4/19/62 


66 


0.40 

Platform 
Displ. 
in 
Inches 


-0.40 


Figure  Xllm  -  Platform  Dj 

.splacements  for  Wave  D 

Theory  

Experiment  ooo 

Run  DO  Unbraced 

X  *              x    •*  """ 

0.30 


0.40 


Wave  Phase  Angle  at  Front  Legs  in  Degrees 


Run  DOZ  Braced 
■M*  -  0.860 


■0.40 


560 


-0.00 


67 


VCH  4/19/62 


Figure  Xlln  -  Platform  Displacements  for  Wave  D 
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Figure  XIIo  -  Platform  Displacements  for  Wave  D 
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Figure  XIIp  -  Platform  Displacements  for  Wave  D 
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IV.   DISCUSSION  OF  RESULTS 

The  plots  of  Figure  XII  indicate  that  the  results 
of  this  thesis  are  good,  in  the  opinion  of  the  authors. 
The  procedure  set  forth  in  the  thesis  is  Judged  to  be 
sound. 

The  reader  will  notice  that  in  all  cases  the 
theoretical  curve  of  displacement  has  the  same  shape 
as  the  experimental  curve  (which  is  not  actually  faired 
in).   Of  particular  note  are  the  curves  of  runs  AO  and 
AOZ  (Figure  Xlla),  for  which  both  the  theoretical  and 
experimental  curves  have  an  extra  dip  in  the  positive 
portion  of  the  curve. 

An  attempt  was  made  during  the  experiment  to  re- 
late some  reference  time  in  the  movie  and  on  the  San- 
born recorder  paper.  The  procedure  involved  pushing 
a  manual  push-button  as  the  second  hand  of  the  clock 
(one  rotation  ■  10  seconds)  passed  some  known  time. 
Any  human  error  would  appear  as  an  error  in  the  experi- 
mentally determined  phase  angles  0  at  the  front  legs . 
The  resulting  data  showed  that  the  human  error  was  so 
large  that  the  data  had  to  be  discarded.   Consequently 
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the  phase  relationships  of  the  experimental  points  in 
Figure  XII  were  adjusted  to  give  the  best  fit  to  the 
theoretical  curve . 

If  one  follows  a  particular  structure  con- 
figuration (for  example,  6z)  through  its  four  runs 
he  will  see  that  the  greatest  displacements  occur  for 
those  waves  which  give  frequency  ratios  ~  of  about  1.1 
or  smaller.   (Note  that  two  different  displacement  scales 
are  used  for  the  plots . ) 

Figure  XIII  contains  one  point  for  each  run.  Each 
point  represents  the  ratio  of  the  maximum  positive 
theoretical  displacement,  plotted  as  a  function  of  fre- 
quency ratio  ~  for  the  particular  run.  The  figure  is 
intended  to  indicate  the  magnitude  of  the  difference  be- 
tween the  theoretical  and  experimental  thesis  results. 
It  is  seen  that  the  displacement  ratio  is  generally  quite 
different  from  its  ideal  value  of  1.0.  Values  greater 
than  1.0  Indicate  a  conservative  design,  or  overdesign. 
Values  less  than  1.0  indicate  that  the  theory  is  in- 
adequate and  that  the  structure  designed  according  to 
the  theoretical  procedure  would  be  too  weak.   It  is  in- 
ferred in  the  previous  paragraph  that  the  designer  is 
interested  in  waves  for  which  7-7  is  about  1.1  or  less. 
Consequently  the  right-hand  portion  of  Figure  XIII  is 
of  no  further  concern,  as  a  practical  matter.   For  runs 
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Figure  XIII  -  Ratio  of  Maximum  Positive 

Platform  Displacements,  Theoretical/Experimental, 
vs.  the  Frequency  Ratio  g  for  Each  Run  of  the 
Thesis  Results  n 
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in  which  -£j  is  less  than  1.1  multiplication  of  the 
theoretical  displacements  by  1.5  would  insure  that  the 
displacement  ratios  were  always  1.0  or  greater. 

The  theoretical  solution  for  displacement  is  very 
sensitive  to  errors  in  the  structural  natural  fre- 
quency.  Computer  calculations  (not  submitted),  using 
both  the  theoretical  and  experimental  values  of  CO 
obtained  in  this  thesis  (Appendix  C),  indicate  that 
in  some  cases  the  maximum  displacements  so  calculated 
differ  by  a  factor  of  2.  The  computer  results  show 
that,  given  two  frequency  ratios  with  all  other  variables 
the  same,  one  may  not  expect  that  the  frequency  ratio 
nearer  1.0  will  generally  give  the  larger  theoretical 
displacement.   If  a  designer  is  in  doubt  about  the  ex- 
act natural  frequency  of  his  structure,  as  will  be  the 
general  case  in  design  work,  he  should  apply  a  factor 
of  ignorance  of  j5  to  his  calculated  displacements. 
Under  certain  circumstances  this  may  be  extremely  con- 
servative (by  a  factor  of  about  5)>  but  at  other  times 
it  may  be  barely  adequate.   (The  designer  may  wish  to 
apply  an  additional  conventional  factor  of  safety.) 

The  authors  desired  that  the  theoretical  calculations 
should  test  the  method  of  determining  wave  forces  on  the 
structure  and  the  method  for  predicting  the  consequent 
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vibrations.  The  calculation  of  theoretical  <>0  was 

n 

felt  to  be  of  secondary  importance  in  testing  this  thesis 
In  order  to  reduce  the  sources  of  errors  in  the  Results, 
therefore,  the  theoretical  calculations  were  made  using 
the  experimental  values  of  60   (as  determined  in  air). 
Theoretical  values  were  used  for  wave  frequency  co, 
however. 

Additional  discussion  of  results  is  given  in 
Appendices  G  and  I. 
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V.      CONCLUSIONS 

1.  The  procedure  set  forth  in  this  thesis  provides 
a  sound  method  of  calculating  the  dynamic  displace- 
ment of  a  fixed  offshore  structure  acted  upon  by  a 
passing  train  of  regular  waves. 

2.  The  dynamic  displacement  of  the  platform  may  be 
of  interest  per  se  (p.  2).  More  generally,  however, 
the  dynamic  displacement  is  required  for  the  design 
stage  determination  of  the  vibration  stresses  through- 
out the  structure  (pp .  4l -42 ) . 

5.    Field  test  data  on  Texas  Towers  No.  2  (pp.  5-6) 
and  No.  4  (p.  l)  shows  that  the  highest  wave  encountered 
does  not  necessarily  produce  the  greatest  platform  dis- 
placement.  Consequently,  for  design  purposes  one  might 
not  design  an  adequate  structure  simply  by  considering 
only  the  largest  wave  to  be  encountered,  as  has  been  the 
general  past  practice.  One  must  consider  several  waves 
of  various  lengths  (or  periods)  and  steepnesses  over  the 
range  of  possible  water  depths.  The  critical  displace- 
ment-wave combination  must  be  arrived  at  by  an  iterative 
process. 
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4.    For  a  particular  structure  with  its  particular 
natural  frequency  the  greatest  displacements  occur 
for  those  waves  which  give  frequency  ratios  —  of 
about  l.i  or  smaller.  For  these  waves  the  theoretical 
solution  for  displacement  is  generally  conservative, 
with  some  exceptions. 


5.    A  designer  should  apply  a  factor  of  ignorance 
of  3  to  his  calculated  displacement.  Under  certain 
circumstances  this  may  be  extremely  conservative 
(by  a  factor  of  about  5),  but  at  other  times  it  may 
be  barely  adequate.   (The  designer  may  wish  to  apply 
an  additional  conventional  factor  of  safety.) 
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VI.   RECOMMENDATIONS 

L  This  thesis  advances  a  procedure  for  obtaining  a 
"safe"  design  of  a  fixed  offshore  structure.  At 
times  the  procedure  may  lead  to  a  considerably  over- 
designed  (excess  strength)  structure,  which  is  gen- 
erally poor  from  a  cost  viewpoint.   It  is  recommended 
that  future  workers  attempt  to  refine  the  procedure 
in  order  to  reduce  the  amount  of  overdesign  that  may 
be  involved = 

2.   It  is  recommended  that  some  calculations  be  carried 
out  for  some  real  structures  in  the  ocean  and  that  the 
results  be  compared  with  observed  data.  From  this  work 
one  might  confirm  the  validity  of  this  procedure  for 
design  purposes,  and  a  revised  factor  of  ignorance 
might  be  obtained.   For  example: 

a.  Were  the  results  of  the  design  calculations 
for  Texas  Tower  No0  4  adequate?  Could  one  predict  tower 
collapse  for  the  wave  and  wind  conditions  that  existed 
when  the  tower  did  collapse?  How  do  calculated  displace- 
ments compare  for  some  of  the  observed  data  available 
in  [2]  ■> 
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b.  Perform  a  design  calculation  for  the 
Argus  Island  structure  and  compare  the  results  with 
those  of  [19]  . 

c.  Perform  a  design  calculation  for  the 
Buzzards  Bay  Light  Station  structure  and  make  a 
conclusion  regarding  the  adequacy  of  structural 
strength. 

3.  Broaden  the  scope  of  the  design  data  by  extending 
the  limits  of  Figures  III,  IV,  and  V  to  include  soli- 
tary and  breaking  wave  theories. 

4.  Develop  an  analogous  procedure  for  calculating 
platform  translation  and  rotation  for  an  arbitrary 
angle  of  wave  impingement . 

5.  Continue  the  investigation  of  vortex  shedding  on 
a  vertical  cylinder  in  a  wave  train  (Appendix  I).   In 
a  manner  similar  to  that  of  this  thesis  calculate  a 
Fourier  series  representation  of  the  vortex  shedding 
force  and  establish  a  procedure  for  calculating  trans- 
verse platform  displacement. 
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VII .      APPENDIX 
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A.  Details  of  Theoretical  Procedure 
Spring  Constant  k 

Apply  simple  beam  theory  [29J  to  one  leg 
(Figure  XIV)  to  determine  the  static  deflection  curve 
for  a  force  P  acting  at  the  level  of  the  platform. 


^  "  EI  *  El(2  "  Z) 

dy  m   fo/iz   zf v 
HZ    EIV2    2  ' 

y   -  ET(T  "  Z~> 

_PoP 

ymax~  12EI 


5Z2   2Z5 


W  "  K2         V 


Po   :  ,  12EI 


By  definition:  k  =  ^ *  --7-3 


«    «  1       1     12NEI       /.,/rx  2 

For  N  legs:  k  »  -p—     (46)    Flgure  XIV  -  Simple 

*  Beam  Representation 

uBT  One  Leg 


Rayleigh's  Energy  Method  Applied  to  the  Unbraced  Structure 

Apply  Rayleigh's  Energy  Method  ("22,50"]  to  the  unbraced 
fixed  offshore  structure  to  determine  the  undamped  natural 
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frequency  CO    .   Assume  simple  harmonic  motion  for  the 
N-legged  structure.   Use  the  static  deflection  curve 
for  the  elastic  curve. 


tf  2         2 

Maximum  kinetic  energy  -  t^O^  y„J[     +  N   I     |g-a>n  ymax 


"     2  3^ 

3Z  2Z^ 


2 


«2    *5 


dz  -  gg^Bpi  ♦  i|  Nw<?) 


M2  dZ   No  A  ' 

Maximum  potential  energy  ■  N   pgy    "  p4EI 


Equate  the  maximum  kinetic  and  potential  energies  to 
obtain: 
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n   '  J5(W  +  ±2n*J{)  (47) 


Use  (36)  and  (46)  to  find: 


m  -  |(W  +  ^NwJ?)  (41) 


81 


If  the  configuration  had  been  such  that  the  legs 
were  considered  to  be  pinned  to  the  platform  the  fraction 
to  be  used  in  the  second  term  of  (4l)  would  be  different 
from  that  calculated  above.   A  calculation  for  this  new 
value  would  proceed  in  a  manner  exactly  analogous  to 
that  above.   One  can  deduce  that  in  calculating  the 
spring  constant  k  he  would  start  with: 


d2y 


y-|f-tf<*-«> 


The  calculation  would  give  a  fraction  of  y$/\k§ 
to  be  used  in  the  second  term  of  (4l ) . 

Still  Water  Added  Mass  for  the  Unbraced  Structure 

Assume  that  the  added  mass  for  any  elemental 
length  of  leg  is  equal  to  the  mass  of  the  water  dis- 
placed by  that  leg  |~3l")>  measured  in  terms  of  weight 
w.j  per  unit  length.   Calculate  the  maximum  kinetic  energy 
K.E.  of  this  water  for  a  condition  of  free  vibration. 
Then  calculate  how  much  of  this  added  mass  unit  must  be 
placed  at  the  level  of  the  platform  (W  )  to  give  the  same 


maximum  K.E. 
•d 


2  _.  2 
2g 


"w^n  ymax 


o 


5Z2   2Z3)   dZ  =  WwC02  v   2 


J5 


n  •'max 
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Use  this  expression  to  calculate  W  for  d  ■  27  in. 


and  H  *  42  in. : 

Ww  -  0.05  N  Wwjj 


By  comparing  this  quantity  with  those  of  (4l)  one 
sees  that  the  added  mass  effect  is  small.  Added  mass 
will  therefore  be  neglected. 

Still  Water  Viscous  Damping  for  the  Unbraced  Structure 

A  numerical  exanple  must  be  used  for  this  cal- 
culation.  Use  run  0-W,  which  is  the  free  vibration  of 
the  unbraced  structure  in  still  water  with  zero  pounds 


added  to  the  platform  in  the  experiment.  Let/y2  denote 
the  root  mean  square  velocity  of  the  leg  maximum  velocity 
profile.  Assume  simple  harmonic  motion. 

D  -  0.505  in.        d  *  27  in.  J^«  42  in. 


yo2  -  \       \         W|£2    jp 


'./*$-  °-359  *«* 
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This  occurs  at  Z  »  17.0  in.,  which  can  be  found 
when  one  realizes  that  maximum  velocity  at  any  point  Z 
on  the  leg  is  related  to  maximum  displacement  at  that 
point  by  the  constant  o>  •  Enter  Figure  VI  with 
a  -  17.0  in.,  a/Jl    -  0.406,  and  find  Po/P  -  O.360. 

For  run  0-W  the  average  of  the  observed  amplitudes 
for  the  decaying  exponential  displacements  was  0.66  in. 

The  experimental  natural  frequency  was  CO     ■  6. 15  rad/sec 
Assume  q  and  oi)  to  be  numerically  equal,  which  is  very 
nearly  true  at  the  small  values  of  damping  which  are  of 
concern.  With  these  numbers  one  can  compute  an  average 

*max  -  "V  W  -  (6-^)(^#)  -  0.337  ft./sec. 

,-.  /J*  -  (0.359)  (0.337)  -  0.121  ft./sec. 
V  o 

(0.121) (0. 505) 

which  gives  a  Reynold's  number  iR  - ^| —  =  442. 

(1.15  x  10  p) 

Enter  Figure  II  to  find  steady  state    C  m   1.2. 

The  maximum  damping  force  acting  on  the  platform  is 

P  '  2 

_o  N  CD  f  Ajj  yQ      where  A^  -  D.d  and  N  =  4. 

This  maximum  damping  force  is  calculated  to  be 
0.00232  lb.   At  the  low  Reynold's  number  involved  one  can 
assume  that  the  product  of  drag  coefficient  CD  and  velocity 
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is  constant,  which  means  that  the  damping  force  varies 
directly  as  the  first  power  of  velocity  rather  than  the 
second  power.  This  assumption  means  that  the  maximum 
damping  force  equals C  ly^  so  that 

C  -  Js  N(CD.yo)  {  D-^f "  0-0192. 

From  (41)  find   m  -  0.272.   From  (57)  find 

c 

C     ■  2m co     »  3.34.      One  can  calculate  viscous  £—  ■  0.0057, 
en  \j  _ 

c 

which  happens  to  be  exactly  the  value  that  was  obtained 
experimentally.   It  is  11$  of  the  observed  structural  damp- 
ing (Table  X). 

With  expressions  for  C  and  C  one  can  write: 

c 

^o  N(const. )&D.d. (CO  .y  ) 
C     P_ n  °, 

<L  2m  CO 

c  n 


For  the  experimental  runs  all  variables  in  the  right 

hand  side  of  this  expression  are  fixed  except  y  and  m. 

(60   cancels  out.)  If  a  situation  is  considered  in  which 
x  n 

y  is  held  constant,  and  m  is  varied,  it  is  seen  that 

u  o 

n 

viscous  7r-   should  decrease  as  m  is  increased.  This  trend 

°c 
was  observed  experimentally. 
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Influence  Fraction  =— 

Suppose  that  the  force  P  acting  as  shown 
(Figure  XV)  on  a  two-legged  structure  causes  the  plat- 
form deflection  y   .  The  problem  is  to  find  what  Prt 

max  o 

acting  at  the  platform  level  will  cause  the  same  ymax< 


M, 


max 


zzzzzzzzzzzzzzzzzzzzzzzzzza 


T 

a 

1 


ir>niii,\,i  1 1 ,  ) )  1 1 1 1 1  ,  n  i  ,  >>  i  n ,  n 


J^. 


tf  x 


Figure  XV  -  Portal  Representation  of  Unbraced  Structure 


From  (46)  it  can  be  shown  that 


ymax   24EI 


and 


Prt 
P  .  _° 
r2   2 
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From  Figure  XV  the  reader  can  deduce  that: 

PQ 
?1   *  P  -  ?p        (sum  of  forces  «  0) 

and  Pa  =  M1  +  M2  +  -|-       (sum  of  moments  ■  0) 


For  0  ^=-   Z  ^  a 


M(Z)  -  Ft.  -  P-jZ 


!   M.Z2   PnZ5 

'  -  ir  4 H 

For  a  ^  Z  —  i  M(Z)  -  M.  -  P-jZ  .+  P(Z-a) 

,i„L   Mn  "  p-iz  +  P(Z-a) 
y  "  EI    x    -1 

1  P1Z2      72 

y'  s  El   Miz  "  TT"  +  P(l~  '  aZ)   +  C3 

Pa2 
Equate  expressions  for y1  at  Z  -  a  to  find  C,  -  -ngr 


x    ILZ-   PZ'     z3   ^2     a2z 

y'ET   "2 5"""  +P(5 2-)+-2—  +C4 


-Pa5 


Equate  expressions  for  y  at  Z  *  a  to  find  C^  ■  -^= 

P 

Let  y'  -  0  for  .Z-=«jpand  use  p  =  p  -  _° 
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P  0 
to  find  Mx  -  Pa(l  -  |j)  -  -<p- 


Substitute  these  values  for  P-,  and  M1  into  the  ex 
pression  for  y  and  evaluate  ym  .  at  Z  -J?  to  find: 


max 


ymax  *  EI 


Pa2J?   Pa^   PoP 
"%  6     2T 


Equate  this  to  the  previous  value  for  y    from 
•p  max 

(46)  and  solve  for  =°-: 


Po   5a2   2a^ 
P    -p  -  j3 

This  expression  is  plotted  in  Figure  VI 


Fourier  Series  Representation  of  Force  Acting 
on  Platform 


The  following  two  expressions  have  been  developed 
theoretically  to  represent  the  force  acting  on  the  plat- 
form: 

F(t)  -  F  sin2cot  +  F  coscot   (0°^o)t  ^l80°)       (55) 
F(t)  -  -F  sin2 cat  +F3coswt  (l80°^  tot  £  360°)     (56) 

Many  textbooks,  such  as  1*22,301,  discuss  Fourier 
series. 
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Let  P(t)  -  a.-,    sincot  +  a^sin2cut   +  ...    +  a  sin  n«Jt 
l  2  n 

+brt  +  b.coscot   +  b_cos2o)t  +  ...   +  b  cos  ncot 
o         1  2  n 


where:        a     ■  — 
n 


CO 


(0°/  cot   £  360°) 
2w/co 


F(t)     sin  ncot  J     dt 


bn  =  = 


b     ■  ■=— 


2r/i 


•jp    /  F(t)[cos  nuit       dt 


One  can  deduce  that  0  =  ap  ■  a^  ■  a,-  =   ... 

and  0  =  b,   ■  b_,  ■  bc  »   .... 

15         5 


coF. 


w/. 


JO 


a 


1  w 


sin5  cot  dt 


COF, 


IT 


2f/ 


CO 


sin5 out  dt  -  fe(F1+  F2) 


*/ 


OpF. 


r/c 


CO 


*5 


sin  cot  I 


sin5CL)t  dt  - 


60 


OOF. 


2ir/c 


CO 


rr/o 


[aln2cot\  sin?  cot] 


tjTT/60 


-  T5F(pi  +  p2> 
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a5" 

10^(pl  +  p2) 

a7  = 

l59iF<Fi  +  F2> 

a9" 

s§5f<pi  +  p2) 

«/c 


00 


2f/o> 


0>F 


2  ir 


i      1         rsin2cotjrcos26Utldt   -  -^ 


J 


I  sin  cot  |jcos2<otjdt 
*7co 


-(*!   -  F2) 


b4  =  0  -  b6  -  b8  =   


*/{ 


CO 


P21T/C 


fc>F. 


b     - 


o  2ir 


Bin  cot  dt  - 


GOF, 
~2t 


CO 


F      -  P 
2  IP 

sln^cot  dt  -  -^—r, — - 


0    ir/o 


CO 


Fl   -  F2       4 


'  •    F(t)   =     ±  4  +  ^^   +  |   )sincot 


-   ^^(F-l   +  F2)sin5cOt    +  F,  cos  cot 


(pi  -  p2> 


cos2Wt   +   .  . . . 


(57) 
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B.   SAMPLE  THEORETICAL  CALCULATIONS 

The  calculations  performed  to  determine  theoretical 
displacements  caused  by  the  B  wave  are  given  below. 

Given: 

pile  diameter  D  -  0.5  in.  (0.0417  ft.) 
water  depth  d  -  27  in.  (2.25  ft.) 
wave  height  H  »  0.313  ft. 
wave  length  L  «  10.86  ft. 

therefore    d/L  »  0.207  and  H/L  =  O.O289 

Sample  Calculation: 

1.   Hydrodynamics 

a. )  Calculate  maximum  drag  forces  and  centroids  on 
a  single  leg. 

0*  £  9  £  l80° 


drag  force  -  FD  =  £  |  d  CDA   where  ^«  62.4  lbs/ft5 

(12) 


FD  =  6.58  CDA 


from  computer  program  results  (Appendix  D)  or  by  in- 
terpolation of  Figure  IV,   A  -  0.00421 
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for  water  temperature  =  6l°P,   0  «  1.21  x  10~^  ft.  /sec 


07  o   H   fl  +  >    sinh  (2J)  "1   /r 


J-  2r|=  1.30 


slnh   (2J)  =  6.69       ,       slnh  J  -  1.698 

So  p  ML[  1  +  3(3-14)  (0.0289^(6.69)  1  .  Q 
d         2(2.25)L  8    (1.698r  J 

^7o 
1  +  — J=—  also  could  be  read  directly  from  the  com- 
puter program  results  (Appendix  D), 

1  +£2°.  (§1)     »  1.0715 
d    vd  'max        ^ 


/R  .   0-0^17    J  32.2(2.25) (0.00421)  =  i8q4 
1.21  x  10"5  1  1'0715 


from  Figure  II   for  /R  -  1804,  CD  =  O.96 
therefore   FD  -  6.58(0.96) (0.00421)  =  0.0266  lbs 

drag  force  centroid  =  s1 
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From  computer  program  results  (Appendix  D)  or  by 

S ' 
interpolation  of  Figure  III,  -r-  -  0.685 

S'  =  0.685(2.25)  -  1.5^2  ft. 

180°  ^  9  -  56O0 

drag  force  -  (F  )  »  "  ^  D  d   CD  trough  B       (15 ) 


Vt  -'6'58  CD  trough  B 
from  Figure  V  (or  Appendix  D),    B  -  0.00276 


■•— -*Vrfe 


rrjo  ~  H  _mo     h  (51) 

d  d  d  v^    ; 


/?yo  "  H 


0.0715  -  I7IP  -  -0.0677 


or  from  Appendix  Dy         1   +   (-2° )  =  0.9325 


0.0417         7  32.2(2.25) (0.00276)  m  15gi 
/Rtrough       1  ^  x  10-5f  0.9323 


from  Figure  II  for    IRtrough  -  1591,  CD  trough  =0.925 
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therefore  (FD)t  -  -6.58(0.925) (0.00276)  -  -0.01702  lbs 
drag  force  centroid  ■  (§"•  ), 


from  Figure  III  (or  Appendix  D),  *   't  -  0.559 


(s')t  -  0.559(2.25)  -  1.258  ft. 

b.)     Calculate  maximum  inertia  force  and  centroid  on  a 
single   leg. 

Fx-      ^fHtanhJ  (25) 

for  J  -  1.^0        ,       tanh  J  -  0.862 
p     m  62.4(0.0417)2    (5.14)  (0.315)  (0.862)  _  0t022^ 


S 


y  D2  LH         /    1  +  J  sinh  J  -  cosh  J  "1 
I  "  3  fl  cosh  J  u  J 


sinh  J  *  I.698        ,        cosh  J  -  1.971 


S      .  62.4(0.0417)2(10.86)(0.513)  f  1  +  1.3(1.698)   -  1.971? 
I  8    (0.02295) (1.971)  L  J 


S     m  1.26  ft. 
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therefore  (FD)t  *  -6.58(0.925) (0.00276)  -  -0.01702  lbs 
drag  force  centroid  «  (§"').»- 


from  Figure  III  (or  Appendix  D),  '   't  -  0.559 


(sf)t  -  0.559(2.25)  -  1.258  ft. 


b . )  Calculate  maximum  inertia  force  and  centroid  on  a 
single  leg. 

p  -  *€-w  E   tanh  J  (25) 

for  J  -  l.^O   ,   tanh  J  »  0.862 
p  m   62.4(0.0417)2  (5.14)  (0.313)  (0.862)  m  QM^ 

o"    ]L  p2  LH    /  1  +  J  sinh  J  -  cosh  J  "1 
^1  =  8  FI  cosh  J  u  J 

sinh  J  -  I.698   ,   cosh  J  -  1.971 


S   »  62. 4(0. 0417)^(10. 86) (0. 313)  f  1  +  1.3(1.698)  -  1.971 7 
I      8  (0.02295) (1.971)      L  J 


S  -  1.26  ft. 
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c.)  Calculate  wave  period. 


T  -  £  where  C  J  &  *«"»»  Jf1+^f[  2(cosh(2J))2+2coSh(2J)^1? 
c  '         J  (       L'H  8    (sinh  3V  ) 


(8) 


.  _i/32.2(2.25)(0.862)J"l+(?.l4(0.0289))2r2(6.77)2+2(6.77)-^' 
C"K  V^O^  7  Z.      s(1.698)4 ^~J 

C  =  6.98  ft. /sec. 

Therefore   T  -  g*  |  -  1.558  sec. 


6i>  =  |E  .  4.05  rad./sec 


C  is  determined  internally  in  the  computer  program 
(Appendix  D)  but  is  not  one  of  the  outputs.  A  simple 
addition  to  the  "WRITE  OUTPUT,"  "DIMENSION,"  and  "FORMAT" 
statements  would  include  C  as  an  output. 

Summary  of  steps  a,b,  and  c: 


FD  »  0.0266  lbs.      S«  -  1.542  ft. 


(FD)t  -  -0.01702      (S')t  -  1.258  ft. 
295  lbs.     3L  -  1 

u)  -  4.03  rad./sec. 


F-j.  -  0.02295  lbs.     S_  -  1.26  ft. 
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Enter  these  figures  in  the  upper  left-hand  corner 
of  the  tabular  calculation  sheet,  Figure  XVI.   (Actually 
in  this  sample  calculation  the  numbers  do  not  agree  ex- 
actly with  those  of  Figure  XVI,  which  were  determined  at 
another  time,  but  the  slight  differences  are  negligible 
for  purposes  of  illustration. ) 

2.  Vibrations 

Consider  run  B6.  This  is  the  unbraced  model 
with  6  lb.  of  added  weight  on  the  platform,  in  wave  B. 

The  reader  should  trace  the  steps  on  Figure  XVI 
as  he  reads  the  following  explanatory  notes. 

n 

Values  of  co  ,  =— ,  and  k  are  taken  from 
n  cc 
Appendix  F.  Otherwise  they  would  be  determined  as  ex- 
plained in  the  Procedure,  Section  II. 

P 
Values  of  the  influence  fraction  p-^  are  read 

from  Figure  VI. 

The  values  of  0  at  the  front  legs  are  simply 
reference  values.  Values  of  platform  longitudinal 
displacement  are  calculated  for  these  twelve  reference 
values  of  0,  and  a  curve  is  faired  through  the  points 
(as  was  done  in  the  preparation  of  Figure  XII ). 

Below  the  values  of  9  the  corresponding  re- 
ference angles  for  each  of  the  significant  terms  in 
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FIGURE  XTI  -  TABULAR  CALCULATION  FOR  PLATFORM  DISPLACEMENT 


the  Fourier  series  (57)  are  listed.  It  is  apparent  that 
a  cosine  varying  force  term  leads  the  sine  varying  force 
term  of  the  same  frequency  by  a  90°  time  phase  angle. 
In  order  to  fit  the  cosine  terms  to  the  classical  dif- 
ferential equation  of  forces  (32)  90°  is  added  to  each 
of  the  cosine  term  reference  angles,  after  which  the 
term  is  treated  as  a  sine  function. 

In  the  middle  portion  of  the  table  the  actual 
values  of  the  angles  for  each  of  the  terms  are  listed. 
For  the  front  legs  the  angle  0   is  determined  from  (40) 

or  Table  II  for  each  value  of  ^S^f   where  n  is  the  co- 

^n 
efficient  of  60 1.  The  angle  0   is  subtracted  from  the 

reference  angle  for  each  case,  because  displacement 

lags  force  (33).  At  the  rear  legs  the  force  lags  the 

force  at  the  front  legs  by  the  lag  angle  n.jZLr  (58 ) . 

For  each  specific  case  the  angle  at  the  rear  legs  is 

found  by  subtracting  h#RL  from  the  corresponding  angle 

at  the  front  legs . 

In  the  lower  portion  of  the  table  the  values  of 

x  (33)  are  calculated  for  each  term  of  the  Fourier 

series  for  displacement.  The  expression  for  the  force 

coefficients  F  are  given  in  (57)=  Values  of  the  mag- 

x 
nification  factor —  are  determined  from  (38) 


XSTATIC 


or  Table  I 
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The  reader  is  cautioned  that  linear  interpolation 

of  Tables  I  or  II  may  be  poor  in  the  immediate  vicinity 

of  resonance  (-rr-  *  1.0). 

^n 

In  each  of  the  twelve  columns  the  value  of  x  is 

o 

multiplied  by  the  sine  of  the  appropriate  angle  from 
above,  and  the  result  is  recorded.  These  results  are 
added  in  each  column  to  give  values  of  platform  longi- 
tudinal displacement  x. 

For  this  sample  calculation  the  ratio  of  the  maxi- 
mum positive  theoretical  displacement  to  the  maximum 
positive  experimental  displacement  is  1.17. 
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C  Summary  of  Theoretical  Calculations 

Forces  and  Force  Centrolds 

Table  IV  describes  the  4  waves  used  in  this  thesis. 
The  values  H  and  L  are  measured,  while  — -=■  '°   is  read 
directly  from  the  computer  output  in  Appendix  D  (or  from 

(6 ) ) .   For  all  runs  d  ■  2 .  25  f t .   Frequency  (radians/ 
second)  is  calculated  from  (8)  by  use  of  the  relation 

6U=  2irC/L.-. 

TABLE  IV  -  Wave  Characteristics 

Wave  H(ft.)  L(ft.)   d/L   H/L  J=2rd/L  ^ffi  d+l|0"  1    co 

A  0.240  15.25  0.1476  0.01575  0.928  1.0557  0.949  3-10 

B  0.313  10.86  0.207  O.0289  1.300  1.0715  0.932  4.00 

C  O.361  8.82  0.255  0.0409  1.603  1.0817  0.921  4.60 

D  0.404  7.33  0.307  0.0552  1.930  1.0908  0.911  5.13 

Table  V  lists  parameters  used  for  computing  drag 
forces.  The  multipliers  A  and  B  are  read  directly  from 
the  computer  output  (or  from  Figures  IV  and  V,  re- 
spectively), IR  and  trough  are  computed  from  (18)  and 
(19),  respectively,  while  Cp  and  CD  trough  are  read  from 
Figure  II. 
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TABLE  V  -  Various  Coefficients  for  Drag  Force  Determination 


Wave 

A 

B 

IR 

CD 

'^trough 

Q 

D  trough 

A 

0.00272 

0.00191 

1410 

0.846 

1188 

0.846 

B 

0.00421 

0 . 00276 

1804 

0.925 

1591 

0.885 

C 

o . 00558 

0.00518 

2070 

0.925 

1720 

0.925 

D 

0 . 00674 

0.00^42 

2290 

0.962 

1795 

0.925 

Table  VI  lists  the  theoretical  forces  and  force 
centroids  on  a  single  leg  of  the  structure.  FD  and  (FD). 

are  computed  from  (12)  and  (15),  respectively,  while  Fj 
is  computed  from  (25).  The  centroids  S'  and  (S '  ).   are 
read  directly  from  the  computer  output  (or  from  Figure  III), 
and  Sj.  is  computed  from  (28).  Centroidal  distances  are 
measured  from  the  bottom  of  the  tank. 


TABLE  VI  -  Theoretical  Forces  and  Force  Centroids  on  a 


Single  Leg 

- 

Wave 

FD(lbs.) 

S'(ft.) 

(FD)t(lbs.) 

(3T')$(ft) 

FJ(lbs) 

Sl(ft.) 

A 

0.01755 

1.575 

-0.01242 

1.180 

0.01485 

1.200 

B 

0.0266 

1.542 

-0.01702 

1.258 

0.02295 

1.260 

C 

0.0557 

1.680 

-0.0205 

1.550 

0.0285 

1.528 

D 

0.0425 

1.815 

-0.0216 

1.408 

0.0529 

1.580 
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Structure  Natural  Frequencies 

Theoretical  natural  frequencies  of  the  structure 
are  computed  by  (47)  and  (6l),  using  numbers  from 
Appendix  E.   In  Table  VII  run  5,  for  example,  designates 
the  unbraced  model  with  5  pounds  of  added  weight  on  the, 
platform.  The  letter  Z  designates  the  braced  model. 
Experimental  values  of  a)   (Appendix  F)  are  also  given 
in  Table  VII  for  convenience  of  comparison.  The  units 
of  all  frequencies  are  radians/second. 

TABLE  VII  -  Structure  Natural  Frequencies 


Run 


0 

5;  25 

5 

4.55 

6 

4.05 

10 

5.59 

oz 

5.20 

5Z 

4.50 

6z 

4.02 

10Z 

5.57 

:p.  in  air 

^n  exp.  in  water 

6.01 

6.15 

4.67 

4.79 

4.00 

4.02 

5.21 

5.27 

5.97 

6.09 

4.95 

5.00 

4.10 

4.11 

5.25 

5.44 
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D.   Computer  Programs 

Two  separate  IBM  7090  FORTRAN  computer  programs  were 
written  for  use  with  this  thesis.   Part  1,  "MAX  DRAG 
FORCE  MULTIPLIER  AND  DRAG  CENTROID,"  was  used  both  to 
develop  the  curves  families  of  Figures  III,  IV,  and  V, 
and  to  supply  exact  values  for  force  and  centroid  de- 
termination for  the  4  waves  studied. 

Part  2,  "TOWER  DISPLACEMENTS,"  was  used  to  supply 
the  theoretical  platform  displacements  for  the  32  ex- 
perimental runs  of  this  thesis.  The  output  from  this 
program  is  shown  plotted  as  Figures  Xlla-p,  The  output 
from  Part  1,  although  comprising  data  for  76  waves  of 
differing  characteristics,  is  included  in  this  Appendix 
only  for  those  waves  which  were  the  object  of  experi- 
mental study  in  this  thesis.   (The  remaining  output  data 
has  been  submitted  to  the  Thesis  Supervisor  under 
separate  cover,) 

Because  the  designer  may  want  to  make  use  of  the 
programs,  a  brief  description  of  their  characteristics 
is  included. 

Part  1 

This  program  was  developed  to  facilitate  solutions 
of  drag  forces  and  drag  force  centroids  for  crest  and 
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trough  regions.   It  is  applicable  only  for  those  wave 
parameters  for  which  Stokes8  J>rd   approximation  applies 
Its  primary  use  is  for  development  of  curve  families 
such  that  A,  B,  -*—  and  v  ,J        in  the  formulas  listed 
below  may  be  read  directly  by  the  designer  if  no  com- 
puter is  available., 

Pd  _  r  ™2z  ci»  a  (12, 

-  y  D  d2  C^  . 
(FD)t  =  — g  D  troU*h  B  (15) 


1+  %/d         2 

S»       /^     Si"  amax  d(S'/d) 

d~  *  d    gd 


A 

2 
u. 


B 


1+Yo/d 

where  As     /     umax.d(S'/d) 

gd 


(26) 


(S8)t       /     S±       trough  d(S8/d)  ,p7n 

d   =  d  ~~ga  ^n 


(14) 
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1-f^o-H 
d 

.2 


Utrough.d(S!/d)  ^j 


gd 


Integration  of  (26,27,14,  and  15)  Is  performed 
by  the  trapezoidal  rule  with  21  unequally  special 
stations. 

The  curve  families  developed  (Figures  III,  IV, 
V)  show  lines  of  constant  H/L  on  graphs  whose  co- 
ordinates are  the  desired  values  vs.  d/L. 

A  secondary  function  of  the  program  is  to  supply 
horizontal -component  velocity  profiles  at  wave  angles 
9  =  90°  and  270° „  Two  uses  of  these  profiles  are  to 
give  specific  indication  of  the  limit  of  applicability 
of  the  authors'  design  procedure  and  to  supply  required 
information  for  an  analysis  of  transverse  vortex 
shedding  forces  (Appendix  I). 

Supplementary  information  which  may  be  read  directly 
from  the  output  includes  s 


d  +  rfo       d  +  Cyp  -  H)   Umax    Utrough 
d    '       d       '  /id   '    /id 
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Wave  celerity  C,  although  determined  internally  as 

— — ,  is  not  included  in  the  output .  However,  the 
v/gd 
simple  inclusion  of  this  parameter  (called  "T"  in  the 

program)  in  the  WRITE  OUTPUT  TAPE,  DIMENSION,  and  out- 
put FORMAT  statements  would  cause  it  to  be  written 
output „ 

Variables  to  be  read  into  the  program  on  data 
cards  are  defined  as  follows: 

X  «  d/ft 
Y  «  L/H 

The  above  data  for  each  run  is  put  on  the  same 
card  in  accordance  with  FORMAT  statement  10.  Variables 
written  as  output  ares 


X  *  d/H 

Y  -  L/H 

A 

CENTC  -  |p 

B 

CENTT  -  (^}t 

PC  ■  (S8/d)n 
UC  -  (-*&)n 

n  »  0,  20 
n  -  0,  20 
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u  2 
UCSQ  -  ("||^)n  n  -  0,  20 

PT  -  (^§Hfih)n         n  -  0,  20 


UT  -  (trough,  n  -  0,  20 

^T  n 

u     2 
utsq  .  (trough^       n  «  0,  20 


Values  of  PC,  UC,  UCSQ,  PT,  UT,  and  UTSQ  are  given 
in  incremental  steps  of  S/d  (not  S'/d)  of  0.05  from 
the  bottom  to  the  still  water  level.   At  n  =  20  it  is 
seen  that  the  following  results 

tL   -  — g 

U 
UC  =  J*** 


ygd 

U  2 
UCSQ  *- 


PT  - 


gd 
d  ±   ("[o   ~   H) 


UT  .  trough 
/gd 

2 

OTaft_  .trough 


The  output  is  written  in  accordance  with  FORMAT  state- 
ments 90  and  100,  and  is  included  in  this  appendix  for 
the  A,  B,  C,  and  D  waves  only. 
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If  other  than  J>6   input  data  cards  are  to  be  used, 
the  program  must  be  run  separately  for  each  set  of 
56,  with  blank  cards  being  used  to  fill  out  any  set 
not  comprising  36  data  cards.  A  card  with  zeros 
punched  in  the  x  columns  should  be  used  to  separate 
the  data  from  the  blank  cards . 

Some  of  the  interesting  aspects  of  the  output  used 
to  develop  Figures  III,  IV,  and  V  are  discussed  in 
Appendix  G. 
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*  Ml 935-1598, FMS, DEBUG, 2, 4, 1000,0 


*  XECJ 

*  LIST 

W.  C.  NOLAN  AND  V.  C.  HONSINGER 

C     WAVE  INDUCED  VIBRATIONS  IN  FIXED  OFFSHORE  STRUCTURES 
C     PART  I,  MAX  DRAG  FORCE  MULTIPLIER  AND  DRAG  CENTROID 
C 

DIMENSION  X(36),Y(56),A(56),B(56),CEN!HC(56),CENTT(36},PC(21,56), 
1PT(21,56),UC(21,3d),UT(21,56),UCSQ(21,5o),UTSQ(21,56) 

READ  INPUT  TAPE  4,10,    (X  (I )  ,Y(l)  ,1-1,36) 
10  FORMAT    (2F10o4) 

DO  60  1-1,36 

N-I 

IF    (X(I))    70.70,20 

20  Z«6.2832*X(I)/Y(I) 
D»EXPF(Z)-EXPF(-Z) 
E*TANHF  (Z ) 

F-EXPF  (2 .  *Z  )-HSXPF  (-2 .  *Z  ) 

T«SQRTF((E/Z)*(l„+(9.8697A(l)**2.)*(F**2.+2.*F+10.)/b**4.)) 
R-0 

DO  50   J-1,21 
GP-EXPF (Z  *R ) +EXPF ( -Z *R ) 
GM«EXPF(Z*R)-EXPF(-Z*R) 
HP-EXPF (2 . *Z *R )  -ffiXPF ( - 2 . *Z *R } 
HM*EXPF (2 0  *Z  *R ) -EXPF ( -2 .  *Z  *R ) 
PP=R+9>4248*HM/(2.*X(I)*Y(I)*D**4. ) 

UU-(9.8697*T/((Y(l)**2.)*D**2. ))*(-2.+6.*HP/D**2.+HP) 
PC (J,I)«PP+GM/(2,*X(I)*D) 
PT(J,l}-PP-GM/(2.X(l)*D) 
UC(J,l)-UU+3.l4l6«T*aP/(Y(l)*D) 
UCSQ(J|l)«UC(J,l)**2o 
UT(J,I)-UU-3.l4l6*T*GP/(Y(l)*D) 
UTSQ(J,I)-UT(J,I)**2. 
IP(J-I)    30,30,40 
A(l)-0 


+  (UCSQ ( J, I )+UCSQ ( J-1,I ) )* (PC ( J, I )-PC (J-1,I )  )/2. 
+  (UTSQ(J,I)+UTSQ(J-l,l))*(PT(J,l)-l 


PT(J-l,l))/2, 
WC  +  (UCSQ(J,I)"+UCSQ(J-1,I))*(PC(J,IJ**2.-PC(J-1,I}**2,  )/4. 
WT+(UTSQ(J,I)+UTSQ(J-l,l))*(PT(J,l)**2.-PT(J-l>l)**2.)/4. 


WC 

WT 
50  R»R+,05 

CENTC(1)-WC/A(I) 
60  CENTT(I)*WT/B(I) 
70  DO  80  M«1,N 

WRITE   OUTPUT  TAPE   2,90,X  (M),Y(M)  ,A(M),CENTC  (M)  ,B(M),CENTT  (M), 
1PC(1.M),UC(1,M),UCSQ(1,M),PT(1,M),UT(1,M),UTSQ(1,M) 

WRITE   OUTPUT  TAPE   2,100, (PC (J,M),UC (J,M),UCSQ(J,M),PT(J,M) ,UT( J,M 
l),UTSQ(J,M),    J-2,21       ) 

80  CONTINUE 

90  FORMAT  (2F10<i6,10F9.6) 
100  FORMAT (56X,6F9„6) 
CALL  EXIT 
END  (1,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0)  109 


I/H 


d/H 

.180000] 63.500000 


0.002722 


0.611706|0. 001910" 


7.180000  1 34. 700000  10.004206  10.6854  31  10.00  27  59  10.559225  |o. 


0.052390 
0.104785 
0.157191 
0.209614 
0.262060 
0.314534 
0.367043 
0.419591 
0.472186 
0.524833 
0.577539 
0.630310 
0.683153 
0.7360  75 
0. 789083 
0.842185 
0.095387 
0.948699 
1.002129 
1.055685 


16.230000  124.450000  10.005381  10.746012  0.003177  10.590742 


5.575000  18. 140000  0.006742  0.806107  0.003418  0.624576 


0.052704 

0.105420 

0.  158160 

0.210935 

0.263759 

0.316644 

0.369602 

0.422646 

0.475791 

0.529050 

0.582438 

0.6359  70 

0.689663 

0.743532 

0.797597 

0.851875 

0.906386 

0.961152 

1.016196 

1.071540 

0. 

0.052724 

0.105465 

0.158243 

0.211074 

0.263979 

0.316974 

0.370081 

0.423321 

0.476713 

0.530282 

0.5840  50 

0.638043 

0.692287 

0.746812 

0.801648 

0.856828 

0.912388 

0.968365 

1.024800 

1.081740 


0.052578 
0. 105180 
0.157831 
0.210556 
0.263380 
0.316330 
0.3694  35 
0.422724 
0.476227 
0.529978 
0.584013 
0.638371 
0.693093 
0.748225 
0.803816 
0.859920 
0.916597 
0.973913 
1.031938 
1.090752 
110 


gd 
0.042524 
0.042578 
0.042740 
0.043011 
0.043391 
0.043882 
0.044486 
0.045204 
0.046040 
0.046994 
0.040072 
0.049277 
0.050612 
0.052082 
0.053692 
0.055449 
0.057358 
0.059425 
0.061659 
0.064067 
0.066658 
0.044294 
0.044401 
0.044724 
0.045263 
0.046024 
0.047009 
0.048225 
0.049680 
0.051382 
0.053342 
0.055572 
0.058087 
0.060901 
0.064033 
0.067505 
0.071337 
0.075557 
0.080192 
0.085275 
0.090841 
0.096930 
0.041733 
0.041884 
0.042338 
0.043100 
0.044176 
0.045575 
0.047310 
0.049396 
0.051853 
0.054702 
0.057970 
0.061688 
0. 065890 
0.070616 
0.075914 
0.081835 
0.088440 
0.095798 
0. 103985 
0.113090 
0.123214 
0.036887 
0.0370/9 
0.037655 
0.038624 
0.039998 
0.041792 
0.044031 
0.046741 
0.049958 
0.053723 
0.058086 
0.063106 
0.068849 
0.075398 
0.082844 
0.091297 
0.100882 
0.111747 
0.  124064 
0.138031 
0.153882 


0.001808 
0.001813 
0.001827 
0.001850 
0.001883 
0.001926 

10.001979 
0.002043 

!0. 002120 
0.002208 
0.002311 
0.002428 
0.002562 
0.002713 
0.002883 
0.003075 
0.003290 
0.003531 
0.003802 
0.004105 
0.004443 
0.001962 
0.001971 
0.002000 
0.002049 
0.002118 
0.002210 
0.002326 
0.002468 
0.002640 
0.002845 
0.003088 
0.003374 
0.003709 
0.004100 
0.004557 
0.005089 
0.005709 
0.006431 
0.007272 
0.008252 
0.009395 
0.001742 
0.001754 
0.001793 
0.001858 
0.001951 
0.002077 
0.002238 
0.002440 
0.002689 
0.002992 
0.003361 
0.003805 
0.004341 
0.004987 
0.005763 
0.006697 
0.007822 
0.009177 
0.010813 
0.012789 
0.015182 
0.001361 
0.001375 
0.001418 
0.001492 
0.001600 
0.001747 
0.001939 
0.002105 
0.002496 
0.002886 
0.003374 
0.003982 
0.004740 
0.005685 
0.006863 
0.008335 
0.010177 
0.012487 
0.015392 
0.019053 
0.023680 


trough :   trou«h 
d  gd 

f-O. 040005 
.047752-0.040039 
.095500hO. 040144 
.  143239H0. 040319 
.190965-0.040563 
.238673-0.040877 
.286359-0.041262 
.334019-0.041717 
.3816481-0.  042242 
.4292411-0.042837 
•  476794J-0. 043503 
.5243031-0.044240 
.571762-0.045047 
.619166-0.045926 
.666512-0.046875 
713795J-0. 047895 
0.048985 
-0.050145 


u 


LZ 


2 

trough 


283846 
33099C 
378057 


518681 
565322 


750  36 
796172 


761008 

808148 

.855209-0.051375 

.902186-0.052675 

949075J-0. 054043 

■0.043083 

04737CJ-0.  043160 

09473C-0. 043393 

14206E-0. 043780 

•l8937«-0. 044323 

236637(-0. 045024 

■0.045882 

-0.046900 

•0.048078 

425034-0. 049419 

471915-0.050925 

-0.052596 

-0.054435 

611824-0.056444 

658175-0.058623 

704 36q-0. 060974 

•0.063497 

0.066192 


84176S-0. 069059 


0.072095 
0.075296 
0.041139 
0473151-0.041254 
094613-0.041598 
14187e-0. 042173 
-0.042981 
236236-0.044024 
283295J-0.  045306 
33025^-0.046830 
3770811-0.  048601 
423769^-0.050623 
470292(-0. 052902 
0.055442 
-0.058249 
0.061327 
0.064680 
-0.068310 
-0.072220 
0.076407 
833786J-0. 080869 
877724-0.085597 
0.090581 
-0.036639 
047440-0.036791 
094856^0. 037246 
0.038008 
0.039081 
0.040470 
2837981-0.042182 
330726|-0.  044227 
•0.046614 
424020J-0.  049354 
470  325-0.052458 
516357-0.055937 
562080-0.059801 
607456-0.064062 
.652442-0.068725 
696994-0.073796 


88713S 
93226 


51663] 
562765 
608666 
654312 
699675 
744729 
78944 


921227 


14222 

189521 

236721 


741063 


.827536 


0.079275 


784596-0.085155 


-0.091421 


869820-0.098046 


0.104989: 0 


z&r 

.001600 
00160  3 
.001612 
.001626 
.001645 
.001671 
001703 
001740 
.001784 
.001835 
.001893 
.001957 
002029 
002109 
.002197 
.002294 
002400 
002515 
002639 
002775 
002921 
001856 
001863 
001883 
.001917 
001965 
002027 
002105 
002200 
002312 
.002442 
.002593 
.002766 
002963 
003186 
003437 
003718 
004032 
004381 
004769 
005198 
005670 
001692 
001702 
001730 
001779 
001847 
001938 
002053 
002193 
002362 
00256- 
002799 
003074 
003393 
003761 
004183 
004666 
005216 
005838 
006540 
007327 
008205 
001342 
001354 
001387 
001445 
001527 
001638 
001779 
001956 
002173 
002436 
002752 
003129 
003576 
004104 
004723 
005446 
006285 
007251 
008358 
009613 
011023 


PART  2 

This  program  performs  the  tabular  calculation  of 
Figure  XVI.   It  is  written  for  the  specific  case  of 
N  *  4  legs  and  k  ■  0.625  lb. /in.   In  the  expressions 

for  S,  T,  U,  V,  and  W  the  constants  are  recognized 

N     N1 
to  be  r  or  r-0   For  different  values  of  N  or  k  the 

constants  may  be  changed  by  substituting  new  values 

in  the  program. 

Variables  to  be  read  into  the  program  on  data 

cards  are  defined  as  follows: 

y  *  run  number 

A  «  co/co 
'      n 

B  -  C/C0 

C  -  0RL   (radians) 

DD  -  -i-rj — -,  (lbs.) 

.  4   (F.  +F2)    (lbs.) 

o  .    i_  (pi  +  p2)     (lbs-) 

M  15w 

R  -  F,    (lbs.) 

The  above  data  for  each  run  is  put  all  on  the 
same  card  in  accordance  with  FORMAT  statement  10. 


XXI 


Variables  written  as  output  are: 

y  »   run  number 

x  m   platform  longitudinal  displacement  (in. ) 

This  output  is  written  on  a  single  line  for  each  run, 
in  accordance  with  FORMAT  statement  100.  For  each  in- 
put data  card  12  values  of  platform  displacement  are 
computed,  these  values  corresponding  to  the  12  angular 
wave  positions  9  indicated  at  the  top  of  Figure  XVI. 
If  more  than  52  runs  are  to  be  used,  the  program 
must  be  run  separately  for  each  set  of  J>2,   with  blank 
cards  being  used  to  fill  out  any  set  not  comprising 
32  data  cards. 


112 


*M1 933-1598, FMS, DEBUG, 2, 4, 500,0 

*  XEQ 

*  LIST 

W .C .NOLAN  AND  V .C .HONSINGER 

C     WAVE  INDUCED  VIBRATIONS  IN  FIXED  OFFSHORE  STRUCTURES 
C     PART  2,  TOWER  DISPLACEMENTS 

DIMENSION  Y(32),A(52),B(32),C(32),DD(52),P(32),Q(52),R(52),THETA(4 
1,12),X(32,12) 

READ  INPUT  TAPE  4,10,  (Y(l ),A(I ) ,B(l ) ,C (I ),DD(l),P (I ) ,Q(l ) ,R(I ) ,1-1 
1,32) 
10  FORMAT  (IF  3-0,7F8.5) 

DO  90  1-1,32 

D-2o*C(l) 

FF-2.*A(l)*B(l)/(l.-(A(l))**2.) 

IF(FF)^20o30,30 
20  FF-FF+3. 14159 
50  CONTINUE 

F-ATANF  (FF) 

GG-4.*A(l)*B(l)/(l.-4.*(A0;))**2.) 

IF(GG)  40,50,50 
40  GG-GG  +  3.14159 
50  CONTINUE 

G-ATANF  (GG) 

HH«6,*A(I)*B(I)/(1.-9.*(A(I))**2.) 

IF(HH)  60,70,70 
60  HH-HH  +  5 . 14159 
70  CONTINUE 

H-ATANF    (HH) 

AA-1./SQRTF((1.-A(I)**2.)**2.    +   (2.*A(l)*B(l ) )**2. ) 

BB-1./SQRTP((1.-4.*(A(I})*»2.)**2.    +     4.*A  I  •Bjl)   **2.) 

CC*l»/SQRTF((lo-9.*(A(I))**2.)**2.    +   (6.*A (I )*B(I ) )**2. ) 

S-6 . 4*DD (I ) 

T«3.2*P(I)*AA 

U-5«2*Q(I)*CC 

V«5o2*RfI)*AA 

W—  3.2*DD(I)*BB 

DO  80  J»l,12 

Z-J 

THETA (1 , J )-0 . 5236 *Z 

THETA(2,J)-1.5708*Z 

THETA  h,  J  )«1. 5708  +  0.5256*Z 

THETA (4, J )-1.5708  +  1.0472*Z 

X(I,J)«S  +  T*SINF(THETA(1,J)-F)  +  U*SINF(THETA(2,J)-H)  +  V*SINF(TH 
1ETA(3,J)-P)  +  W*SINF(THETA(4,j)-G)  +  T*SINF (THETA (1, j)-P-C (I ) )  +U 
2*SINF(THETA(2,J)-H-E)    +  V*SINF(THETA (5, J)-P-C (I ) }   +  W*SINF(THETA(4 

5,J)-Q-D) 

80  CONTINUE 
90  CONTINUE 
100  FORMAT (1F6o0,12F7. 3) 
DO  110  1-1,32 

WRITE  OUTPUT  TAPE  2,10C>Y(I ) ,  (X(I,  J),  J-1,12) 
110  CONTINUE 

CALL  EXIT  ,,-, 

END (1,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0)  **2 


E.   Sample  Calculations  for  Experimental  Data 

B  wave 

Wave  period  T: 

Sanborn  recorder  speed  ■  25  mm/sec. 

Wave  length  as  measured  on  recorder  ■  40  mm, 

Therefore   T  »  40/25  *  1.60  sec. 

Wave  circular  frequency  co  : 

2ir  2ir 


CO 


T"  "  I75o  "  5*95  radlans/sec 


Wave  height  H  is  read  from  the  experimental  wave  profile 
(Figure  XVTIb). 

Wave  forces: 

Wave  drag  and  inertia  forces  are  read  directly  from 
the  experimental  force  record  (Figure  XVIIb ) . 


F    (9  -  90°)  •  0.044  lbs. 


(PD)t    (0  -  270°)  -  -0.039  lbs. 
Fj    (0  -  180°)  -  0.017  lbs. 


Properties  of  the  Structure  Model 

Consider  run  6  and  calculate  the  structural 
damping  by  (42)  and  (44). 
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flVPrflffP     1  -  1/.92    ,    .93       .64       .62       .43\       n   k~ 
average  —  -  ^(-^  +  33  +  77^  +  753  +  73$)  -  1#47 

cf-  loge(l.47)  -  O.385  (42) 

§-  -  ^L  -  0.0613  (44) 

c 

Damped  natural  frequency: 

q  -  4?73ysecS(2ir)  ~  5*99  radians/second 
From  (45): 

a)  •  —    ^'"       «  4.00  radians/second 
1 1   -  (0.0613)2 

The  volume  of  plastic  in  the  unbraced  model  (in- 
cluding base)  is: 

Lai  )2  -  »(7  )2|o.6lo]+  i,5(0.505)2(42)  +  (50^  -  £)(22)(|) 

-  199  +33.6  +  329  -  561.6  in5 
For  one  leg: 

wJl   -  (|)(|ol^)(22|)  4:0.340  lb 


'IB.  lb: 
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For  the  platform: 


The  bolt  weighs  0.200  lb 


,199 


3\ 


W  -  (§fl7o-)(22^)  +,0.200'^)8.25.1b. 


i  j 


Calculate  the  following  non-dimensional  ratios 
From  (59) t 


xm              (ii 

^(0.540) 

From    (63): 

PD  -  1  + 

149 

1 

L  2J 

(4)  (27) 

J 

From    (65 ) : 

rx-l  + 

149 

"        ~1 

1 

1 

i2\ 

(*){ar) 

1.69 


-  1.34 


%l§ 


F.  Summary  of  Calculations  for  Experimental  Data 

Wave  Characteristics 

Wave  profiles  for  waves  A-D  are  shown  in  Figure  XVII 
a-d,  respectively.  The  crest  is  in  all  cases  the  8th 
wave  crest  to  pass  the  cylinder.  Wave  phase  angle  9 
is  set  equal  to  90°  at  the  crest. 

TABLE  VIII  -  Experimental  Wave  Characteristics 


Wave 

H(ft.) 

L(ft.) 

d(ft.) 

T 
(sec . ) 

CO 

(radians/sec ) 

A 

0 .  240 

15.23 

2.25 

2.00 

3.14 

B 

0.315 

10.86 

2.25 

1.60 

3.93 

C 

0.361 

8.82 

2.25 

1.40 

4.49 

D 

0.404 

7.33 

2.25 

1.26 

4.99 

Longitudinal 

Wave  Forces 

Experimental  and  theoretical  longitudinal  forces 
associated  with  waves  A-D  are  shown  in  Figure  XVIIa-d, 
respectively.  The  forces  are  read  as  follows: 


FD  at  9  «  90° 


(FD)t  at  9  -  270' 

Pj  at  9  -  0° 


Figure  XVIIa  also  shows  experimental  transverse 
forces  caused  by  vortex  shedding  in  the  A  wave. 

The  experimental  longitudinal  forces  read  from 
Figure  XVII  are  given  in  Table  IX. 


1" 
TABLE  IX  -  Experimental  Wave  Forces,  ^  Cylinder 


Wave 

FD  (lbs.) 

(Vt  (lbs.) 

FI  (lbs.) 

A 

0.014 

-0.019 

0.017 

B 

0.044 

-0.039 

0.023 

C 

0.060 

-0.055 

0.022 

D 

0.085 

-0.039 

0.032 
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Figure  XVIIa 
Longitudinal  and  Transverse  Forces 
l/2n  Cylinder  -  A  Wave 
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Figure  XVI lb 
Longitudinal  Forces 
1/2"  Cylinder  -  B  Wave 


o  o  o  Experimental 


Wave  Profile  -  B  Wave 
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Figure  XVI I c 
Longitudinal  Forces 
1/2"  Cylinder  -  C  Wave 
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Wave  Profile  -  C  Wave 
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Figure  XVTId 
LONGITUDINAL  FORCES 
1/2 *  Cylinder  -  D  Wave 
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Properties  of  the  Structure  Model 


TABLE 

X  -  Experimental  Damping 

and  Structure 

Natural  Frequencies 

Run 

Total 

Structural 

Viscous 

Damped 

Undamped 

Damping 

damping 

damping 

Natural 

Natural 

(as  fraction 

C 

(as  fraction 

Frequency 

Frequency 

of  critical) 

*c 

of  critical) 

q (rad/sec ) 

uJ  (rad/sec ) 

0 

0.0501 

0.0501 

0 

6.00 

6.01 

3 

0.0548 

0 . 0548 

0 

4.66 

4.67 

6 

0.0613 

O.O613 

0 

3.99 

4.00 

10 

0 .  0748 

0 . 0748 

0 

3.20 

3.21 

0-W 

0 . 0558 

0.0501 

0.0057 

6.12 

6.13 

3-W 

0.0580 

0.0548 

0.0032 

4.78 

4.79 

6-W 

O.O655 

O.0613 

0.0042 

4.01 

4.02 

10-W 

0.0748 

0 . 0748 

0 

3.26 

3.27 

OZ 

0.0501 

0.0501 

0 

5.96 

5.97 

3Z 

0.0559 

0.0559 

0 

4.94 

4.95 

6Z 

O.0613 

O.0613 

0 

4.09 

4.10 

10Z 

0 . 0738 

0.0738 

0 

3.24 

3.25 

oz-w 

0.0547 

0.0501 

0.0046 

6.08 

6.09 

3Z-W 

0.0592 

0.0559 

0.0033 

4.99 

5.00 

6Z-W 

0 . 0645 

O.0613 

0.0032 

4.10 

4.11 

10Z-W 

0.0730 

0 . 0738 

■" _ 

3.43 

3.44 

Spring  constant  k  (46)  was  determined  experimentally: 
k:  =»  0.625  pounds/inch 
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G.  Supplementary  Dlscusaion  of  Results 
1.  Theory  and  Experiment 

Longitudinal  Forces 

A  glance  at  the  longitudinal  force  plots  of 
Figures  XVII  a-d  shows  considerable  disagreement  be- 
tween theory  and  experiment .  Theoretical  forces  (drag 
and  inertia)  are  tabulated  in  Table  VI,  and  experimental 
forces  are  in  Table  IX.  These  tables  show  close  agree- 
ment of  inertia  forces,  with  the  drag  forces  differing 
by  a  factor  of  2,  in  cases. 

Prior  users  of  the  force  gage  have  indicated  that 
the  dynamic  response  of  the  gage  to  exciting  forces 
may  indicate  forces  greater  than  the  actual  /_28J. 
However,  the  authors  conclude  that  for  the  frequencies 
associated  with  the  waves  used  in  this  thesis  any  magni- 
fication factor  applied  to  force  gage  displacements 
would  in  all  cases  be  less  than  1.1. 

The  proposal  to  use  the  -pms  velocity  in  the  Reynolds 
numbers  for  drag  coefficient  determination  is  based  upon 
the  assumption  that  C^  and  C-  ^ToURh   are  independent  of 
depth  j17j.  For  the  waves  studied  the  lowest  velocity 
encountered  in  the  crest  region  is  (theoretically)  at  the 
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bottom  of  the  tank  for  the  D  wave.   From  Appendix  D 

it  is  seen  that  max  ■  O.0369  for  this  situation. 

Therefore,       's    u    -  O.315  ft. /sec.  and  /R  = 
'  max     •-'/'         » 

1092  for  this  region.  When  this  value  is  compared  with 
the  Reynolds  numbers  listed  in  Table  V,  a  glance  at 
Figure  II  shows  that  the  above  assumption  is  valid  for 
the  waves  studied  in  this  thesis.  However,  if  a  local 
|R  at  the  bottom  of  a  structure  leg  becomes  less  than 
about  100  in  ocean  waves,  the  assumption  will  lead  to 
inaccuracies. 

The  drag  and  inertia  components  were  of  the  same 
order  of  magnitude  for  the  experimental  runs  of  this 
thesis.  Harleman  and  Shapiro  J"  27I  have  found  that  for 
this  particular  situation  agreement  is  not  good  between 
the  proposed  theoretical  forces  and  the  experimental 
forces.  Because  the  theoretical  and  experimental  inertia 
forces  of  this  thesis  agree  quite  closely,  (Tables  VI 
and  IX),  it  is  suspected  that  the  discrepency  is  due  to 
the  method  of  evaluating  the  drag  force  component.  Ex- 
periments conducted  by  Wiegel,  Beebe,  and  Moonr5lto 
determine  average  values  of  CD  in  a  wave  train  show 
much  scatter,  with  CD  varying  between  0.2  and  2.5  for 
Reynolds  numbers  of  around  10  .  They  suggest  that  a 
factor  of  between  2.0  and  2.5  be  applied  to  the  drag 
force  if  an  average  value  of  CD  is  used. 
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Keulegan  and  Carpenter [yzl ,   in  their  work  with 

horizontal  cylinders,  have  related  C-.  to  the  period 

u   T 
parameter  max   rather  than  to  some  Reynolds  number. 

D 

If  the  drag  coefficient  is  to  be  averaged  over  the 

height  of  a  vertical  cylinder,  the  authors  suggest  that 
CD  and  (CT))trouKh  be  included  with  A  and  B  in  the  com- 
puter program  (Appendix  D).  A  look-up  table  for  drag 

u   T 
coefficient,  based  upon  the  curve  of  Cn  vs  max 

v  D 

in  \j>2 J,  would  have  to  be  included  in  the  program. 


Free  Oscillations  and  Natural  Frequencies 
of  the  Structure  Model 


In  Appendix  H  it  will  be  seen  that  the  free  os- 
cillation data  is  somewhat  odd  in  that  when  the  structure 
swings  from  a  negative  peak  to  a  positive  peak  the  ampli- 
tude of  displacement  often  increases.  Theoretically  for 
the  damped  system  this  amplitude  must  decrease [22] . 
When  the  authors  saw  this  result  they  immediately  sus- 
pected unequal  spring  constants  in  the  two  directions, 
but  a  test  showed  that  the  spring  constants  were  the 
same.  The  authors  are  unable  to  explain  this  oddity. 
It  is  seen  in  Appendix  F  (or  C)  that  in  all  con- 
figurations the  natural  frequency  of  the  model  is  higher 
in  still  water  than  in  air.   One  expects  just  the  reverse, 
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because  of  viscous  damping  by  the  water  and  the  added 
mass  effect.  The  authors  cannot  explain  this  result, 
but  they  conclude  that  the  added  mass  effect  is 
negligible . 

It  is  also  seen  that  for  all  runs  with  the  added 
weights  the  frequency  of  the  braced  structure  model  is 
higher  than  that  for  the  unbraced  structure  model,  both 
in  air  and  in  water.  Theoretically  spring  constant  k 
was  not  changed  when  the  model  was  braced.  Consequently, 
the  addition  of  the  bracing  mass  should  have  lowered  the 
natural  frequency  u)  (  -  V  k/m ) .   Because  the  frequency 
is  higher,  however,  for  the  braced  model,  it  is  probable 
that  the  bracing  increased  the  stiffness  k  slightly.  The 
authors  did  not  suspect  this  result  at  the  time  of  the 
experiment,  however,  and  consequently  stiffness  k  was 
not  measured  for  the  braced  model. 

2.   Computer 

It  is  mentioned  in  Appendix  D  that  velocity  profiles 
at  crest  and  trough  regions  were  written  as  output  for 
the  computer  program  used  to  develop  Figures  III,  IV  and 
V.  These  profiles,  most  of  which  are  not  submitted  in 
this  thesis  (see  comment,  p.  103),  show  some  interesting 
aspects  relative  to  the  boundary  limits  on  Figure  V.  The 


127 


profiles  of  particular  interest  are  those  for  the 
trough  region.  The  computer  output  gives  values  of 

u         S ' 

trough      trough 

i&  d 

For  certain  critical  values  of  wave  parameters 
d/L  and  H/L,  the  peak  value  of  rtrough  is  shown  to 
occur  at  some  level  other  than  at  the  surface,  some- 
times at  the  ocean  bottom.   In  the  extreme  situation 

utrough  at  the  surface  ("  utrough^  becomes  positive. 
Table  XI  gives  the  values  of  ut rough  associated  with 

these  critical  parameters.   The'^  dimensionless 

velocity  trough  refers  to  the  velocity  at  the  surface, 


and 


V  gd      bottom  refers  to  the  velocity 
at  the  ocean  bottom    '  s    under  the  wave  trough. 

O  I 

The  dimensionless  parameter  (  trough )     refers  to 

d    '  max   S '     . 
the  surface  elevation  of  the  trough,  and   ( r°ug  )cv±+ 

refers  to  the  elevation  at  which  the  absolute  velocity 
is  either  a  maximum  or  zero,  as  appropriate.  Values  of 
d/H  and  L/H  are  included  in  Table  XI  because  the  com- 
puter output  is  identified  by  these  terms. 

The  authors  suggest  that  the  reader  refer  to 
Figure  V  while  interpreting  the  information  in  Table  XI. 

First,  for  d/L  -0.07  and  H/L  -  0.04,  it  is  seen 
from  Table  XI  that  the  predicted  height  of  the  trough  is 
above  the  still  water  level.  This  result,  then,  can  be 
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TABLE  XI  -  Trough  Velocities  Associated  with  Critical 


Wave  Parameters 

d/L 

H/L 

d/H 

L/H 

trough 
>/gd 

S  ' 

/     trough \ 

*       d          'max 

oottom    t 
Vsd 

1 

trough \ 
d        'crit 

0.20 

0.10 

2.000 

10.00 

-0.1809 

0.777 

-O.I636 

0.680 

0.25 

0.14 

1.786 

7.14 

-0.1937 

0.739 

-0.1592 

0.650 

0.10 

0.05 

5.333 

33.33 

-O.O96O 

-- 

-0.1057 

— 

0.15 

0.08 

1.875 

12.50 

-0.1588 

-- 

-0.1904 

— 

0.07 

0.04 

1.750 

25.00 

-K3.1269 

1.030 

40.0021 

1.030 

0.08 

0.05 

1.600 

20.00 

-to .  0788 

0.976 

-0.0646 

0.650 

0.10 

0.08 

1.250 

12o50 

+0 . 0685 

0.901 

-0.1686 

0.775 

discarded  as  being  outside  the  range  of  possibility.   It 
is  interesting  to  note  that  this  situation  is  reached  when 
Stokes1  3rd  approximation  predicts  positive  velocities  in 
the  trough  region  from  surface  to  bottom. 

If  the  remainder  of  the  data  is  compared  with  the 
curve  families  of  Figure  V,  the  following  observations 
are  noted: 

1.   Each  curve  of  the  curve  families  reaches  a  peak 

value,  which  occurs  approximately  where  Stokes' 

3rd  approximation  predicts  a  positive  velocity 

S ' 
for  trough  regions  above  some  ( *joug  )crif 
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2.   Values  of  the  parameters  corresponding  to  the 
above  situation  are  seen  to  fall  approximately 
on  line  H/d  «■  O.653,  which  has  been  stated  to 
be  the  theoretical  limit  for  wave  characteristics 
in  the  ocean.   The  profile  of  positive  velocity, 
beginning  at  the  surface  and  extending  down  as 
d/L  is  decreased,  is  a  result  of  the  increasing 
significance  of  mass  transport  in  Stokes'  j5rd 
approximation . 

5.   Values  of  the  parameters  for  which  u.    .  is 
never  positive  (and  for  which  its  peak  value 
does  not  occur  at  the  surface)  are  seen  to  fall 
approximately  on  the  line  for  which  H/d  ■  O.605. 
This  line  corresponds  to  the  limit  of  applica- 
bility of  Stokes'  theory. 

With  regard  to  the  statements  made  above,  it  is 
to  be  noted  that,  when  plotted  on  the  curve  families 
of  Figures  III,  IV,  and  V,  lines  of  constant  H/d  converge 
for  low  values  of  d/L.   Also,  in  explanation  of  ob- 
servation number  1  above,  it  must  be  remembered  that  in 
solving  for  the  force  multiplier  B  from  (15)  it  is 
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assumed  that  u,    ,  is  negative  for  all 
trougn      ° 


trough 
~~cl 


If  the  proposed  procedure  is  to  be  extended  to  regions 


for  which  this  is  not  true,  then 

2 


u 


should  be  substituted  for  u 


trough 


trough 
in  (15). 


u 


trough 


151 


H.  Original  Experimental  Data 

Water  temperatures  (°F)  for  the  various  runs 
are  given  in  Table  XII » 

TABLE  XII »  Water  Temperatures  (°F) 


Wave 

Unbraced 
Structure 

Braced 
Structure 

Forces  on  , 
Single  Cylinder 

A 

62 

57 

65 

B 

62 

61 

63 

C 

62 

61 

63 

D 

62 

61 

6? 

A  typical  wave  gage  calibration  curve  is  shown  in 
Figure  XVIII.  A  typical  force  gage  calibration  curve 
is  shown  in  Figure  XIX.   "Record  reading"  on  these  curves 
refers  to  the  deflection  of  the  stylus  needle  on  the  San- 
born recorder.  The  wave  gage  is  calibrated  about  the 
still  water  level,  while  the  force  gage  is  calibrated 
about  the  zero  force  position.   Typical  Sanborn  recorder 
wave  profile  traces  are  shown  in  Figures  XX  and  XXI. 

The  damped  natural  frequency  of  the  0.5  in.  dia. 
cylinder/force -gage  combination  was  4.8  cps,  while  for 
the  1  in.  dia.  cylinder/force -gage  combination  (used 
for  transverse  forces)  it  was  17.85  cps. 
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Figure  XVIII 
Typical  Wave  Gage  Calibration 
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Amplifier  Attenuation  X20 
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WCN  4/13/62 


Figure  XIX 
TYPICAL  FORCE  GAG 


0/ 


Longitudinal  Force  Calibration 
for  1/2"  Cylinder 


Amplifier  Attenuation  X10 


WCtf  4/13/62 
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Wave 
Profile 


Long. 
Force 


Figure  XX  -  Sample  Sanborn  Record,  Wsve  Profile 
And  Longitudinal  Forces 


Wave 

Profile^ 


Trans • 
Force 


Figure  XXI  -  Sample  Sanborn  Record,  Wave  Profile 
And  Transverse  Forces 
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Properties  of  the  Structure  Model 

Experimental  data  for  free  (natural)  oscillations 
of  the  model  is  given  in  Table  XIII.  Run  3,  for  example, 
designates  a  run  in  air  with  the  unbraced  model  with  3 
pounds  of  added  weight  on  the  platform.  The  letter  Z 
designates  the  braced  model.  The  letter  W  indicates  runs 
with  the  model  in  still  water.  The  tabulated  numbers  are 
peak  displacements  as  obtained  from  movies.  A  number  may 
not  represent  the  actual  peak,  which  might  have  occurred 
between  frame  exposures .  After  viewing  the  movies  the 
authors  deduced  that  in  those  few  cases  in  which  an  error 
could  exist  it  could  be  only  about  0.01  inch.  The  elapsed 
time  listed  is  between  the  first  and  last  recorded  peak- 
displacements. 

The  weight  of  plastic  in  the  unbraced  model  (in- 
cluding base)  is  22-2/4  lb.  The  weight  of  plastic  bracing 
is  about  0.4  lb.  The  weight  of  the  bolt  is  0.200  lb. 

The  total  length  of  bracing  is  186  inches .  The 
length  of  bracing  X-q   below  the  still  water  level  is 
149  inches. 
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TABLE  XLII-   Experimental  Data  for  Free  (Natural) 
Oscillations  of  the  Structure  Model 


Run: 

0 

3 

6 

10 

0-W 

3-W 

6-W 

10-W 

time: 

3»l4 

4.04 

4.73 

5.88 

3.08 

3.94 

4.70 

5.79 

Peak 

-.89 

-.90 

-.92 

-.81 

-.98 

-.96 

-.90 

-  .89 

Displ. 

.98 

.99 

.93 

.83 

.95 

.91 

.89 

•^87 

-.68 

-.66 

-.64 

-.51 

-.70 

-.68 

-.60 

-.56 

.70 

.68 

.62 

•  51 

.67 

.62 

.59 

.54 

-.48 

-.46 

-.43 

-.31 

-.49 

-.48 

-.40 

-.35 

.51 

.49 

.43 

.32 

.47 

.43 

.40 

.34 

-.56 

-.32 

-.29 

-.20 

-.35 

-.33 

-.26 

-.19 

Run: 

oz 

3z 

6Z 

10Z 

oz-w 

3Z-W 

6Z-W 

10Z-W 

Time: 

5.16 

3.82 

4.61 

5.81 

3.10 

3.78 

4.60 

5.49 

Peak 

-.67 

-.69 

-.84 

-.83 

-.97 

-.94 

-.88 

-.78 

Displ. 

.74 

.67 

.98 

.92 

.71 

.93 

.87 

.74 

-.48 

-.48 

-.58 

-.50 

-.67 

-.65 

-.59 

-.49 

.54 

.47 

.66 

.57 

.50 

.63 

.58 

.51 

-.56 

-.34 

-.39 

-.32 

-.47 

-.46 

-.39 

-.30 

.41 

.34 

.46 

.37 

.36 

.43 

.39 

.33 

-.25 

.24 

-.26 

-.18 

-.34 

-.32 

-.26 

-.18 
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I.  Vortex  Shedding  on  a  Vertical  Cylinder  In  a  Wave  Train 

General 

It  has  been  mentioned  in  this  thesis  that 
in  the  development  of  a  practical  offshore  structure  de- 
sign transverse  or  "lift"  forces  caused  by  vortex  shedding 
may  be  included.  It  is  to  be  expected  that  these  reactive 
forces  are  periodic  in  nature,  with  magnitudes  and  fre- 
quencies being  functionally  related  to  the  horizontal  com- 
ponent of  wave  orbital  particle  velocity  u.  Furthermore, 
since  u  *  0  at  discrete  values  of  wave  phase  angle  9,  it 
may  also  be  expected  that  this  periodicity  in  vortex 
shedding  Is  repetitive  with  the  passage  of  each  wave  in 
a  regular  train. 

As  an  illustration  let 


F.  .  ■  magnitude  of  transverse  force  caused  by  a  vortex 
J   being  shed 

9.  .  ■  wave  phase  angle  corresponding  to  F,  . . 

where   i  *  1,2,.. „n  ■  sequential  number  of  a  wave 

in  a  regular  wave  train  (1  =•  1 
for  the  first  wave,  etc.) 

j  ■  l,2,...m  ■  sequential  number  of  the  vortex 
shed  during  the  passage  of  one 
wave  (j  -  1  for  the  first  vortex 
shed  during  the  passage  of  each 
wave,  etc.) 
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Then 

F   ■  constant  (1  ■  l,n) 
■*•  J 

0.  .  ■  constant  (i  ■  l,n) 
but  since  P..  Is  a  function  of  Uj., 
P.  .  ^   constant  (j  «  l,m) 

Thus  it  is  seen  that,  because  of  the  repetitive 
nature  of  these  forces,  if  their  magnitudes  and  fre- 
quencies can  be  predicted  a  Fourier  series  representation 
could  be  utilized  in  developing  a  procedure  for  calcu- 
lating transverse  platform  displacements.  This  procedure 
would  not  differ  appreciably  from  that  presented  in  this 
thesis  for  longitudinal  wave  forces.  It  is  the  intent  of 
this  appendix  to  present  an  approach  for  the  determination 
of  these  vortex  forces  and  shedding  frequencies  and  to 
show  the  results  of  experiments  performed  by  the  authors 
to  substantiate  this  approach. 

A.  Theory 

Transverse  Forces 

Before  the  transverse  forces  are  discussed, 
it  might  be  well  first  to  make  some  qualitative  remarks 
concerning  the  significance  of  considering  these  forces 
in  a  structure  design.   Mention  has  been  made  (Introduction) 


139 


of  the  reported  failure,  due  to  fatigue,  of  a  piling  sub- 
jected to  transverse  vibrations  in  ocean  waves. 

It  is  to  be  noted  also  that  the  frequency  of  this 
force  will  be  at  least  equal  to  wave  frequency,  and  for 
very  long  ocean  waves  it  could  be  many  times  greater 
than  wave  frequency.  It  is  possible  that  this  force 
could  create  a  resonance  situation  in  an  offshore 
structure . 

Nothing  has  been  found  in  the  literature  relative 
to  the  magnitudes  of  these  forces  as  a  vertical  cylinder 
in  a  wave  train,  and,  to  the  authors'  knowledge,  no 
measurements  of  the  forces  have  been  made  in  surface 
waves.  However,  if  one  relates  the  magnitude  of  this 
force  to  the  maximum  drag  force  by  applying  other  type 
measurements  and  steady  state  theory,  some  interesting 
observations  can  be  made. 

Let 

FL  -  maximum  transverse  force  acting 
on  vertical  cylinder  of  length  X 

F-PJ  ■  maximum  longitudinal  drag  force  acting 
on  vertical  cylinder  of  length  J( 

Then 

2 


dF, 


CT  %  D  Umax  d^  (66) 


L  2 


dFD  -  CD  <|  D  umax  d^  (67) 
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In  steady  state  flow  Landweber  [55]  deduced  that 

CL  -  3.44  ^>   4  CD,  while  for  a  vertical  cylinder  moving 

longitudinally  with  simple  harmonic  motion  in  a  stationary 

liquid  Laird  [^41  found  CL  «  0.37  by  measurement.  Laird's 

measurement  was  for  a  cylinder  diameter  of  1"  and  with 

umax  "  lo7  ft«/s'ec.    Laird,  Johnson,  and  Walker  [  J>$1 

found  experimentally  that  (L   0.3  CD  for  a  2"  diameter 

cylinder  moving  under  the  same  circumstances  as  mentioned 

above.  If  C^  is  taken  as  approximately  equal  to  1,  this 

observation  substantiates  the  value  of  CL  »  O.37.  The 

above  values  of  CL  represent  the  extremes  which  the 

literature  gives  as  bases  for  comparing  transverse  and 

drag  forces. 

It  is  suspected  that  the  value  of  CL  given  by 

Laird  most  closely  approximates  the  real  situation. 

However,  two  important  factors  must  be  considered.  First, 

the  artifice  employed  by  Laird  in  approximating  a  vertical 

cylinder  in  surface  waves  is  deficient  in  that  the  water 

velocity  past  his  cylinder  is  not  a  function  of  the  distance 

from  the  water  surface.  In  the  real  situation  not  only  is 

umax  Continually  varying  from  bottom  to  surface,  but  also, 

1 
as  a  wave  crest  passes  a  cylinder,  the  length  over  which 

any  force  acts  is  not  the  water  depth  d  but  rather  d  +°?o 

C^o  -  height  of  wave  crest  above  the  still  water  level). 
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Secondly,  the  experimental  value  of  CL  «  O.37  was  for  a 

given  cylinder  diameter  D  and  given  umax«   It  is  to  be 

expected  that  CL,  like  CD,  is  a  variable  with  dependence 

on  some  defined  Reynolds  number. 

Prom  (66)   it  is  seen  that  dFT  ■  f(u  f).  From  the 
x  L    x  max' 

definition  of  u  „v  it  is  seen  then  that  dFT  is  that  in- 
max  L 

cremental  force  which  would  result  if  a  vortex  were  shed 
at  the  instant  that  the  crest  passed  the  cylinder.  Gen- 
erally the  vortex  is  not  shed  exactly  under  the  crest. 
However,  if  one  is  interested  in  the  maximum  force  which 
might  occur,  as  is  the  designer,  he  can  justify  an  as- 
sumption that  the  vortex  is  shed  under  the  crest,  where 
the  particle  horizontal  velocities  are  maximum.  There- 
fore, it  is  assumed  that  a  vortex  is  shed  when  0  «  90°. 

Because  waves  of  finite  height  are  being  considered, 
it  must  be  noted  further  that  the  design  force  for  the 
crest  region  will  not  be  equal  in  magnitude  to  the  design 
force  for  the  trough  region. 

As  a  further  complication  it  was  noted  by  the  authors 
in  the  experimental  phase  of  this  thesis  that  some 
vortices  shed  in  the  crest  region  were  swept  back  over  the 
cylinder  by  velocity  reversal  during  trough  passage.  These 
cause  significant  additional  forces  on  the  cylinder,  but 
how  they  add,  detract,  or  interfere  with  new  vortices  about 
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to  be  shed  is  a  matter  of  conjecture.  By  the  time  these 
vortices  arrive  back  at  the  cylinder,  some  of  their  vorticity 
has  been  diffused  out  of  the  vortex  street.  Thus  it  is 
difficult  to  determine  if  the  reduction  in  force  in  the 
trough  region  is  the  result  of  new  vortices  being  shed 
in  a  region  of  lower  orbital  velocity  or  if  the  re- 
duction is  caused  by  vorticity  diffusion  of  the  swept- 
back  vortices.  Visual  observation  must  be  made  at. the 
same  instant  that  force  is  recorded  to  resolve  this  issue. 

In  the  ocean,  where  for  long  waves  u  is  in  the 
same  direction  for  a  relatively  long  period  of  time,  it 
is  proposed  that  most  forces  resulting  in  the  wave  trough 
will  be  caused  by  new  vortices  being  shed,  with  only  re- 
latively few  crest-shed  vortices  being  swept  back  by 
velocity  reversal  -  these  few  being  immediately  after 
9  =  l80°.  The  circumstances  imposed  by  laboratory 
techniques  where  short  wave  periods  are  used  make  it  im- 
practical to  verify  predications  on  transverse  forces 
in  way  of  wave  troughs .  The  authors  feel  that  field  tests 
must  be  employed  for  this  type  of  verification. 

For  this  reason,  it  is  suggested  that  if  a  value 
for  FT  for  the  wave  crest  region  is  predicted  and 
verified  experimentally,  then  a  value  for  the  trough  region 
equal  in  magnitude  to  KFr  (K  *-  l)  could  be  used  in  the 
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design  stage.  It  is  proposed  further  that  for  the 
structure  displacement  solution  a  forcing  function 
of  the  following  type  be  used. 

0°  S  9  ^  180° 

P  -  FL  sintO  t  (68) 

180°  £:  9  ^  360° 

P«KFL  sin  ">t  t  (69) 


where    ^  ■  vortex  shedding  circular  frequency  in 
crest  region. 

u).  -  vortex  shedding  circular  frequency  in 
trough  region. 

K  -  constant  (K  <  l)  whose  value  is  dependent 
upon  the  assumed  origin  of  the  force. 


More  will  be  said  about  shedding  frequencies  later. 
The  sine  function  results  from  the  assumption  [_22J  that 
a  vortex  shedding  force  is  varying  sinusoidally. 

The  remainder  of  thifs  section  on  forces  will  be 
concerned  with  the  prediction  of  F,  and  K,  where  K  will 
be  that  value  based  upon  the  assumption  of  new  vortices 
being  shed  in  the  trough  region.  Several  references  may 
be  cited  (for  example  [j>J>, >4,36, j57j  which  can  be  used  to 
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derive  values  of  CL  by  the  steady  state  Von  Karman 
vortex  street  analysis,  and  this  procedure  will  not  be 
repeated  here.  Rather,  the  previously  specified  values 
of  CL  -  3.44  and  CL  -  O.37  will  be  used  to  place  limits 
on  the  predicted  force  in  formula  (66).  C,  ■  5.44  must 

Li 

be  considered  to  be  very  conservative,  however,  since 

its  value  is  based  upon  the  assumption  that  u  =  u  _ 

max 

for  0°  £   9  <l80°. 

From  (66 ) ,   FL  s/bj,  \  D  uj[  d X  (70 ) 

where  the  limits  of  integration  are  as  yet  un- 
specified. Before  attempting  to  solve  this  equation 
for  FT ,  it  might  be  well  to  first  consider  what  criterion 

Li 

may  be  used  to  indicate  the  boundaries  for  which  vortex 
shedding  forces  are  to  be  considered.  For  steady  flow 
a  regular  Karman  vortex  street  in  the  wake  of  cylinders 
has  been  observed  only  in  the  range  of  Reynolds  numbers 
|r(«  ^)   from  about  60  to  5000.  At  lower  values  of  /R 
the  flow  is  laminar,  and  no  vortices  are  shed.  For 
higher  |R  there  is  complete  turbulent  mixing  [l8j. 
For  the  non-steady  flow  associated  with  wave 
passage,  Keulegan  and  Carpenter  [_J2j  have  developed  a 
more  realistic  dimensionless  parameter  called  the  period 
parameter.  From  experiments  with  horizontal  cylinders 
located  at  various  depths  in  a  sinusoidal  wave  train, 
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they  found  that  the  lower  boundary  for  vortex  street 
development  was  coincident  with  a  value  of  the 
"period  parameter"  between  12.5  and  15. 


where  period  parameter  -  B££  (71) 


u   T 

nar 


u    «  horizontal  particle  velocity  under 
crest  at  the  cylinder  depth 

T  *  wave  period 


'   The  significance  of  the  period  parameter  lies  in 
the  critical  distance  a  particle  must  travel  in  one 
direction  for  a  given  velocity  in  order  to  produce  a 
vortex.  For  a  complete  discussion  of  this  approach 

the  reader  is  referred  to  the  original  work  fj>2j. 

umaxT 
In  this  thesis  the  conservative  value  of  • — ~ —  ■  12.5 

will  be  used  for  the  lower  boundary  of  vortex  shedding, 
while  it  is  noted  that  Keulegan  and  Carpenter  observed 
regular  shedding  for   g^  as  high  as  110. 

Two  alternatives  are  therefore  proposed  for  supplying 
the  limits  of  integration  to  (70).  First,  it  may  be  con- 
sidered that  the  vortex  extends  from  the  surface  down  only 

u   T 
to  the  depth  where  —fr*  =  12. 5 -     Secondly,  it  may  be  con- 
sidered that  the  vortex,  once  started  on  the  surface, 

extends  all  the  way  to  the  bottom  of  the  cylinder.  This 

u   T 
second  proposal  presupposes  that  myc  ■*  12.5  is  the 

critical  value  for  starting  a  vortex  on  the  surface,  but, 
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once  the  starting  inertia  is  overcome,  the  period  para- 
meter is  no  longer  critical  in  cutting  off  the  vortex 

at  a  specified  point  along  the  cylinder. 

U   T 
Regardless  of  which  proposal  is  followed,  mgx  -  12.5 

(U  QY  »  uw,QV  at  surface)  affords  the  limiting  boundary 
max    max 

below  which  no  vortices  are  shed.  Intuitively  the 
authors  feel  that  the  true  solution  is  a  compromise  be- 
tween the  two  proposals  outlined  above.  However,  both 
proposals  will  be  tried  for  the  specific  waves  of  this 
thesis . 

For  descriptive  purposes  the  solutions  proposed 
thus  far  are  designated  as  follows: 


(a)  vortex  extends  to  cylinder  bottom. 

d  +°?o 
PT  -    C  CT  %  d  u  2  d^  fQr  crest  region   (72) 


f 


L  L  2    max 

y 

d-(H-7o) 

FT  -  CT  r  D  u^    .  d  J    for  trough 

L  L  2    trough     region  "*  (75) 


where  y  is  the  distance  of  the  bottom  of  the  cylinder 
above  the  bottom  of  the  tank  (or  ocean). 


14? 


u   T 
(b)  vortex  extends  only  to  a  depth  where   mf*  =12.5 

u     ,  T  D 

and  ^-2 —  =12.5  for  crest  and  trough  regions,  respectively 


d+/7°     o      2 
FL  =     C  cl  2  D  UmaX  d/^  for  crest  region     (74) 

y 

d-(H-fo) 
FL  -     C  CL"2  D  utrough  d   for  trouSn  ^gion   (75) 

y 

u   T 
where  y'  is  that  depth  where  mgx  =12.5  and 

u     ,  T 
trough  „  125  for  ^j  and  ^5^  respectively. 


(1)  CL  -  3-44   (Landweber  -  steady  flow)  (76) 

(2)  CL  »  O.37   (Laird  -  simple  harmonic  motion)       (77) 


As  an  example, 

solution  (a-2)  refers  to  (72)  and  (73),  using  a  value 
of  CT  from  (77).   As  mentioned  previously,  it  is  expected 
that  the  true  solution  lies  between  (a-2)  and  (b-2)  with 
solutions  (a-l)  and(b-l)  being  the  conservative  extremes. 

The  solutions  of  (72,73,74,75)  are  of  the  following 

form: 

FL  =  CL|u2i  (78) 
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where 

u  is  the  mean  square  velocity  over  the  length  to 
be  considered.   If  the  vortex  and  the  cylinder  extend 
to  the  bottom  of  the  water  u  can  be  determined  directly 
from  (20)  and  (21). 


-2   .  A  g  d2  (  0) 

umax   d  +09Q  v 

2 
utrough  *  d  +  {%Q   -  H;  (21' 

where  A  is  determined  from  Figure  IV  and  B,  from  Figure  V", 


*- 


rn       3  ir       sinh 


Tm(1+*m    (sinhV     f  <6> 


The  parameter  -jr-  may  also  be  read  from  the  computer  out- 
put (not  plotted),  Appendix  D.  If  the  vortex  does  not 
extend  to  the  bottom  the  following  simplification  may  be 

used.  Reference  fill   shows  that  for  distances  removed 

2 
from  the  bottom  u  vs.  S  is  approximately  linear,  where  S 

is  the  distance  above  the  bottom.  A  conservative  value 

_2 
of  u  may  then  be  taken  as 

-2  .  U2  +U2(UT/D-  12.5)  (7g) 
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This  equation  requires  use  of  (10)  and  (11)  or  the  com- 

_2 
puter  output  in  Appendix  D  only  two  times  to  get  u  . 

The  authors  feel  that  any  loss  of  accuracy  is  justified 

when  the  procedure  as  a  whole  is  viewed. 

The  proposed  solutions  therefore  reduce  to  the 

following  forms: 

(a-1)  crest:   FL  -  3.44  |  D  u^jTd  +  ^Q  -  yj       (80) 

trough:   FL  =  5.44  I  D  u£pough[  d  -  (H  -JQ)   -  yj   (8l) 

(a-2)  crest:   PL  *  0.37  {  D  u^  [&  +  ftfQ   -  y  J      (82) 

trough:   FL  -  0.37  f  D  ^voa&k[&   -  (H  -  %)  -  yj   (83) 


Z^tt2 


(b-l)  crest:   FL  »  ^{0^4  +  ^  -  y'j 


(84) 


trough:   F1  «  3.44  £  D  u^^  d  -  (H  -  njQ)   -  y «}   (85) 


(b-2)  crest:   FL  »  0.37  f  D  u^axf d  +  %  "  y'7 


(86) 


trough; 


:   FL  =  0.37^  D  u^^d  -  (H  -  J0)  -  y'}    (87) 


y  =  distance  of  bottom  of  cylinder  above  water 
bottom 

uT 
y1  ■  distance  from  bottom  to  point  where  ^-  ■  12.5 

and  the  u  •  s  refer  to  the  mean  square  velocity  over  the 

vortex  length. 
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Thei  force  centroids  aLre  exactly  the  same  as  the  drag 
force  centroids  (Figure  III)  if  the  vortex  extends  to  the 

bottom, and  for  any  other  vortex  length  the  centroids  may 

_o 
be  assumed  in  a  manner  analogous  to  the  assumed  u    and 

max 
r2 

L  trough ' 


Shedding  Frequencies 
Let  S  «  cylinder  Strouhal  number  for  steady  flow 

S  -■%£  (88) 

c    2w   u  v   ' 

where  o) 0  is  the  vortex  shedding  circular  frequency  and 
s 

u  is  the  steady  state  velocity.   For  values  of /R  between 

10  and  10  experiments  by  Roshko  /~^J6 1   show  an.average 

value  of  S  of  about  0.206.  This  is  the  range  of  concern 
c 

in  the  experimental  work  of  this  thesis.   For  other  ranges 

of  Reynolds  numbers,  one  is  referred  to  the  reference 

cited  above. 

Using  S.  -  0.206,  one  obtains  for  u)   , 
c  s 


»B   -  ^TT^  (2r)  (89) 


From  (89)  it  is  seen  that  in  a  wave  train  the 
shedding  frequency  (and  number  of  vortices  shed  per  wave 
cycle)  is  a  function  of  the  orbital  particle  velocity  on 
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the  surface  of  the  water  (U).  This  statement  is  made 

because  vortex  shedding  is  initiated  at  the  surface. 

The  question  next  arises  as  to  what  value  to  assign  to  U. 

Laird  P^J    found  that  for  his  cylinders  moving 
in  simple  harmonic  motion  no  vortices  were  shed  until 
the  velocity  reached  70$  of  the  peak  velocity.  If  this 
result  is  to  be  considered  it  seems  reasonable  that  (89) 
may  be  written  as 

.,,  _  0.206(0.7  Umax  +  Umax)(2rr)   0.206(0.85  Umax)(2ir)   ,onx 
^c  2D  D  KJ   ' 

where   \u'(     mag  +  max'  =  "average  velocity"  at  the 

surface  during  vortex 


shedding  in  the  crest 
region 


and 

U 


^  „  0.206(0.85  "trough) (2r)  (9l) 

for  the  trough  region. 


Keulegan  and  Carpenter  £^2j,  on  the  other  hand,  in  their 
work  on  horizontal  cylinders  in  sinusoidal  waves  propose 
that  the  shedding  frequency  be  related  to  the  average 
velocity  over  a  half  cycle.  If  this  reasoning  is  applied 
to  (89)  the  result  is 

w  ^  0.206(0.5  Umax)(2ir)  (92) 

c         D  for  the  crest 
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and 
^  _  0.206(0.5  "trough)  (2ir)  (95) 

for  the  trough 


The  authors  feel  that  it  would  be  reasonable  to  use 
(90)  and  (91)  for  laboratory  work  and  for  steep  waves 
in  the  ocean,  while  (92)  and  (93)  may  be  used  for  long, 
low  ocean  waves  where  the  period  is  long  enough  to  per- 
mit many  vortices  to  be  shed.  Equations  (90)  and  (91 ) 
are  used  by  the  authors  to  obtain  Shedding  frequency. 
The  designer  may  want  to  modify  this  selection,  however, 
when  it  is  realized  that  the  worst  frequency  to  use  is 
that  frequency  which  most  closely  approaches  the  natural 
frequency  u)     of  the  structure.  In  any  case  the  vortex 

force,  of  this  frequency,  would  be  applied  to  the  legs 

UT 
of  the  structure  only  for  the  time  duration  between  g- 

being equal  to  12.5  on  either  side  of  the  crest  and 

trough.  The  forcing  function  proposed  by  the  authors 

assumes  the  force  to  be  active  for  the  entire  wave 

passage. 

Let  (Tfl )  =  time  required  to  shed  1  vortex  in 

c   the  crest  =  2n__  (9b) 

2»c 


(Ta ) .  ■  time  required  to  shed  1  vortex  in  the 

.   t   trough  -  2t (95) 

2u>t 
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Then  from  (90)  and  (91) 


/m  ^  .         D 2.86  D  tr~y. 

(TS)c-2(0.206)(0.85)umax-Tr^-  <*> 

and 

s  L   u trough 

T 

2  ■  time  for  passage  of  wave  crest  or  wave  trough. 

Let 

n  ■  number  of  vortices  shed  during  crest  passage. 
c 

n.  =  number  of  vortices  shed  during  trough  passage 

Then 

_     ^       T _  T     max  fnc,\ 

nc  "  W^     T772V  (98) 

n     -      ,T    y     .  I  „HrP,ugh  (99) 

nt        2(Ts)t        b.72  D  W' 


Because  n     and  n     must  be  integral  numbers,   it  is  pro- 

C  u 

posed  to  accept  the  next  lower  whole  numbers  in  lieu 
of  the  values  computed  from  (98)  and  (99).  These  inte- 
gral values  of  n  and  n  then  give 

C        \> 


6l)c  »  nc^  (100) 

u)t  -  nt">  (101) 

where  cd   is  the  wave  circular  frequency. 
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The  proposed  theoretical  procedure  is  to  use  (80) 
through  (87)  to  determine  transverse  forces,  and  (98) 
through  (101)  to  determine  the  forcing  frequencies. 

B.  Results,  Theoretical  and  Experimental 

For  the  experimental  runs  the  same  waves  (A,  B, 

C,  and  D)  were  used  that  were  used  for  longitudinal 
force  determination.  Transverse  force  measurements 
were  attempted  by  using  the  0.5  in.  dia.  cylinder 
described  in  Section  II -B.  However,  it  was  found  that 
the  B,  C,  and  D  waves  caused  resonant  vibrations  of  the 
cylinder/force -gage  combination  in  the  transverse 
direction,  with  the  A  wave  being  the  only  wave  which 
gave  usable  results.  For  comparative  purposes  the 
transverse  forces  caused  by  the  A  wave  are  plotted  with 
the  longitudinal  forces  and  the  results  are  shown  in 
Figure  XVIIa. 

As  a  consequence  of  these  results  it  was  decided 
to  center  investigation  on  a  1  in.  dia.  cylinder. 
This  cylinder  cleared  the  bottom  of  the  tank  by 
12-5/8".  Therefore,  in  (80,81,82,83)  y  is  taken  as 
12-5/8"  *  1.052  ft.  By  measurement  the  damped  natural 
frequency  of  the  1  in.  dia.  cylinder  arrangement  was 
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17.85  cps,  as  compared  to  4.8  cps  for  the  0.5  in.  dia. 
cylinder  arrangement. 

Before  the  results  are  presented  it  would  be  in- 
teresting to  predict  whether  or  not  vortex  shedding  is 

expected  for  the  various  runs.   Table  XIV  lists  values 

U   T      U,     T 
of  wave  period,  as  well  as  "jg*   and    pUgn  ,  for 

the  4  waves  of  the  experiment.  Wave  periods  are  de- 
termined from  (8),  while  the  velocities  are  taken 
directly  from  the  computer  output  (Appendix  D).   All 
other  values  in  the  following  analysis  are  taken  from 
Appendix  D,  unless  a  particular  equation  is  specified. 

TABLE  XIV-  Wave  Periods  and  Period  Parameters 


Wave 

T(sec. ) 
2.015 

U   T 
max 

D 

U      T 
trough 

D 

A 

13*75 

-11.13 

B 

1.558 

15.45 

-12.00 

C 

lo350 

17*07 

-12.53 

D 

1.200 

18.90 

-12.90 

uT 
Since  the  critical  value  of  ^-  is  considered 

to  be  12.5,  it  is  doubtful  that  the  trough  region 

can  support  vortex  shedding.   In  fact,  all  forces  in 

the  trough  region  were  observed  to  be  caused  (in  part 

at  least)  by  crest-shed  vortices  being  swept  back  to 
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the  cylinder  during  velocity  reversal.  Hence,  trough- 
shed  vortices  are  not  considered  further,  except  to 
state  that  for  the  runs  performed  by  the  authors, 

Parameters  used  in  the  solutipn  of  (80,  82,  84, 

86)  are  shown  in  tabulated  form  in  Table  XV.     The 

parameter  (umax)cpit  refers  to  the  velocity  for  which 

the  period  parameter  is  equal  to  12.5  and  is  used  as 

the  entering  argument  in  Appendix  D  to  determine  y ' . 

The  parameter  (umo  )  is  the  velocity  under  the  crest 

iiioA  y 

at  the  bottom  of  the  cylinder.     The  parameter    (u     2) 

max  a 

is  for  use  with  (80)  and  (82)  and  is  equal  to 

U   2   u   2 

max  +rj  max  &    The  parameter  (u*m,J)b  is  for  use 

with  (84)  and  (86)  and  is  equal  to 

2       2 
"max  x  (W)critt   velocities  are  in  ft. /sec. 

TABLE  XV*  -  Wave  Parameters  for  Transverse  Force  Determination 

Wave  7o<ft>     Umax  (umax>crit  (umax}y  *'<**)  &m£)a  (^J)b 

A  0.120    0.568  0.518  0.399  2.09  0.241  0.296 

B  0.173   0.825  O.669  0.452  2.00  0.442  O.563 

C  0.200   1.050  0.774  0.462  1.97   '  0.657  0.879 

D  0.227  1.310  0.870  0.453  1.95  0.958  1.334 
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The  transverse  forces  (F,)  as  determined  experi- 
mentally are  shown  plotted  with  the  corresponding  wave 
profiles  in  Figures  XXII  a-d.  Table  XVI  compares 
theoretical  and  experimental  results  for  values  of  FL 
and  n   (number  in  parenthesis  is  integral  value).  The 
value  of  «4>   («  m>,  for  reasons  mentioned)  is  deter- 

C        T/ 

mined  from  (100).  The  theoretical  values  of  nc  are 
determined  from  (98),  and  the  various  F*s  from  the 
appropriate  formulas  (80,  82,  84,  86),  with  inputs 
from  Table  V. 

/  -  1.94  slugs/ft.5      d  -  2.25  ft. 


TABLE  XVT  -  Comparison  of  Theoretical  and  Experimental 
Transverse  Force  Results 


Method 


Wave 


A 
B 
C 
D 


a-1 

a-2 

b-1 

b-2 

O.O883 

0.0095 

0.0230 

0.00248 

0.1680 

0.0181 

0.0654 

0.00704 

0.2550 

0.0274  10.1167 

0.01255 

0.3800 

0.0408 

0.1950 

,0.02100 

Exp. 


n 


Theory  Exp. 


0.014 
0.072 
0.128 
0.088 


2.4(2]  2 
2.7(2)  2 
3.0(2]  2 
5.3(2j  2 


a)/2ir(cps) 


Theory  Exp 


0.993  0.993 
1.283  1,283 
2.220  1.482 


2.500 


I.667 


C.   Discussion  of  Results 

From  Table  XVI  and  Figures  XXII  a-d  it  is  seen 
that  the  assumption  of  vortex  shedding  frequency 
being  equal  to  an  integral  multiple  of  wave  frequency  is 
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valid.  Laird  £54]]  found  experimentally  that  vortex 
shedding  appears  to  occur  in  pairs,  and  this  may  be 
why  the  experimental  values  of  n  for  the  C  and  D  waves 
are  equal  to  2  rather  than  3.  For  all  experimental 
runs  trough  shedding  frequency  c<;  was  equal  to  crest 
shedding  frequency  CO    >   &s  explained  previously. 

The  vortex  forces  in  the  trough  region  were  not 
determined  analytically,  because,  for  the  experimental 
runs  of  this  thesis,  these  forces  are  caused  (in  part) 
by  back-sweeping  crest -shed  vortices  from  which  some 
of  the  vorticity  has  been  diffused  out  of  the  street. 
From  Figure  XVTIa. for  the  0.5  in.  dia.  cylinder,  it  is 
seen  that  the  experimental  maximum  vortex  force  is 
4j5. 4$  of  the  experimental  longitudinal  total  force  and 
90$  of  the  experimental  longitudinal  drag  force.  This 
90$  value  is  considerably  higher  than  the  30$  value 
found  experimentally  by  Laird,  Johnson,  and  Walker  [j53j* 

From  Table  XVI  it  is  seen  that  method  (b-l)  agrees 
most  closely  with  the  experimental  force  results.  Method 

(b-l)  refers  to  (84),  where  CL  is  taken  as  3.44,  and 

u   T 
the  vortex  length  terminates  where  max  ■  12. 5 .  Methods 

"IT— 
(a-2)  and  (b-2),  which  assume  CL  -  0.37,  are  seen  to  give 

forces  which  are  far  too  small.   It  is  seen,  therefore, 

that  a  method  of  solution  between  (a)  and  (b)  together 

with  a  method  between  (l)  and  (2),  would  give  analytical 
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results  most  closely  approximating  the  experimental  re- 
sults of  this  thesis.  The  90$  ratio  of  CL/CD,  determined 
experimentally  for  the  0.5  in.  dia.  cylinder  (Figure  XVTIa), 
may  indicate  a  basis  for  selection  of  CL.  More  experi- 
mental tests  are  required  to  narrow  the  boundaries  for 
the  analytical  forces,  as  shown  in  Table  XVI,  but 
it  is  reassuring  that  the  experimental  forces  lie 
within  the  range  defined  by  the  proposed  theories. 

D.  Conclusions  and  Recommendations 

1.  The  proposed  method  of  determining  a  design  vortex 
shedding  frequency  is  valid.  It  is  suggested  that 
the  designer  pick  a  value  of  n  which  is  on  the 
side  of  the  computed  value  to  cause  LQn    (or  0)  ,  ) 

c        1/ 

to  approach  the  structure  natural  frequency  u)  • 

2.  Until  a  more  accurate  prediction  may  be  made  of 
the  vortex  shedding  forces,  a  safe  design  is  to 
assume  C,  -  C~. 

5.   It  is  recommended  that  field  tests  be  conducted 
with  ocean  waves  of  long  period  to  verify  the 
authors •  procedure  for  the  case  of  trough-shed 
vortices. 
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4.   If  platform  rotation  is  considered,  vortex 
shedding  forces  will  assume  added  importance 
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